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Abstract 

In  this  report  the  initiation  and  subcritical  growth  of  surface  cracks  in  graded  materials 
due  to  sliding  contact  are  considered.  After  a  brief  introduction  the  general  coupled  crack/contact 
problem  for  a  semi-infinite  graded  medium  subjected  to  a  sliding  rigid  stamp  of  arbitrary  profile 
is  formulated.  Solving  the  problem  in  the  absence  of  any  cracks,  the  complete  stress  state  on  the 
surface  of  the  medium  is  evaluated  and  critical  stress  that  would  cause  surface  crack  initiation  is 
identified.  The  coupled  problem  is  then  solved,  stress  intensity  factors  are  calculated  and  some 
results  are  presented. 

1.  Introduction 

Graded  materials,  also  known  as  functionally  graded  materials  (FGMs)  are 
multiphase  composites  with  continuously  varying  volume  fractions  and,  as  a  result, 
thermomechanical  properties.  Used  as  coatings  and  interfacial  zones  they  reduce  the 
residual  and  thermal  stresses  resulting  from  the  material  property  mismatch,  increase  the 
bonding  strength,  improve  surface  properties  and  provide  protection  against  severe 
thermal  and  chemical  environments.  Many  of  the  present  and  potential  applications  of 
FGMs  involve  contact  problems.  These  are  mostly  load  transfer  problems  in  deformable 
solids,  generally  in  the  presence  of  friction  as  in,  for  example,  bearings,  gears,  cams, 
machine  tools  and  abradable  seals  in  gas  turbines.  In  such  applications  the  concept  of 
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material  property  grading  appears  to  be  ideally  suited  to  improve  the  surface  properties 
and  wear-resistance  of  the  components  that  are  in  contact. 

From  the  standpoint  of  failure  mechanics  an  important  aspect  of  contact  problems  is 
the  surface  cracking  which  is  caused  by  friction  forces  and  which  invariably  leads  to 
fretting  fatigue.  In  most  applications  material  property  grading  near  the  surface  is  used  as 
a  substitute  for  homogeneous  ceramic  coatings.  In  both  cases  that  is,  in  both 
homogeneous  and  graded  coatings  the  surface  of  the  medium  consists  of  100%  ceramic 
which  is  generally  a  brittle  solid.  Hence,  the  "maximum  tensile  stress"  criterion  may  be 
used  for  crack  initiation  on  the  surface.  Once  the  crack  is  initiated,  its  subcritical  growth 
under  repeated  loading  by  a  sliding  stamp  is  controlled  by  stress  intensity  factors  at  the 
crack  tip.  The  main  objective  of  this  study  is,  therefore,  the  evaluation  of  peak  tensile 
stresses  on  the  surface  for  the  purpose  of  studying  crack  initiation  and  the  stress  intensity 
factors  for  modeling  the  subcritical  crack  growth.  Specifically,  the  objective  is  the 
examination  of  the  influence  of  friction  coefficient  and  material  nonhomogeneity 
parameters  on  the  peak  surface  stresses  and  stress  intensity  factors.  The  problem  is 
considered  under  the  assumptions  of  plane  strain,  Coulomb  friction  and  linear 
nonhomogeneous  elasticity. 

Studies  in  contact  mechanics  in  elastic  solids  were  originated  by  Hertz  [1],  The 
technical  literature  on  the  subject  is  very  extensive.  A  thorough  description  of  the 
underlying  solid  mechanics  problems  in  homogeneous  materials  may  be  found,  for 
example,  in  Johnson  [2].  Some  sample  solutions  for  frictionless  contact  problems  in  a 
semi-infinite  graded  medium  are  given  in  [3]-[5].  Details  of  the  analysis  of  homogeneous 
substrates  with  FGM  coatings  having  positive  or  negative  curvatures  and  extensive 
results  regarding  the  stress  distribution  under  plane  strain  conditions  and  sliding  contact 
are  discussed  in  Guler  [6]. 
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2.  Formulation 


The  coupled  crack/contact  problem  for  a  nonhomogeneous  half-plane  considered  in 
this  study  is  described  in  Figure  1 .  The  half  plane  is  in  sliding  contact  with  a  rigid  stamp 
of  arbitrary  profile.  The  normal  and  tangential  forces  transferred  by  the  contact  are  P  and 
rjP  respectively  where  rj  is  the  coefficient  of  friction,  and  contact  area  extends  from 
y  =  a  to  y  =  b.  The  half-plane  contains  a  surface  crack  of  length  d  which  is 
perpendicular  to  the  surface.  In  this  report,  we  will  formulate  the  problem  and  reduce  it  to 
a  system  of  singular  integral  equations.  Solving  the  integral  equations  numerically  we 
will  examine  the  effects  of  material  nonhomogeneity  and  friction  on  the  stress  intensity 
factors  and  contact  stresses.  Largely,  for  mathematical  expediency  it  will  be  assumed  that 
the  elastic  parameters  of  the  medium  may  be  approximated  by 


Figure  1:  The  general  description  of  the  crack/contact  problem  in  a  graded  medium. 


ix{x)  =  /xoexp(7x),  k  =  constant,  (la,b) 
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where  fx  is  the  shear  modulus,  7  is  the  nonhomogeneity  parameter,  k  =  3  -  Av  for  plane 
strain  and  k  —  (3  —  u)/(1  +  u)  for  generalized  plane  stress,  v  being  the  Poisson's  ratio. 
By  using  the  Hooke's  law 


a„(x,y)=  K[\{(K  +  l)gx+(3  «)-}, 

(2a) 

(2b) 

.  .  ,  .  r  du  dv\ 

^(*,  »)=/*(*){  Q-y  +  dii- 

(2c) 

The  equilibrium  conditions  cr^j  =  0  can  be  expressed  as, 

/  \^2u  ,  . d2u  dv  >du  \®v  n 

<~K+1^  +  (K~1)W+  +^{K  +  1)T*  +  7(3 "  K)di  =  °' 

/  \d1'v  /  O  d2u  1  1\®V  ,  /  i\®u  n 

(k  +  1)v+(k'1)^+  5^  +  T(K"  ^  7(K“  ^  = 


(3a) 

(3b) 


In  previous  studies  (e.g.  Delale  and  Erdogan  [7])  it  was  shown  that  the  stress  intensity 
factors  in  graded  materials  are  not  significantly  influenced  by  the  variation  in  v.  Thus,  in 
this  study  too,  the  Poisson's  ratio  will  be  assumed  to  be  constant.  Following  boundary 
conditions  must  be  satisfied  in  the  solution  of  the  problem 


<?xx(0,  V )  =  0,  axy{ 0,  y)  =  0,  -  00  <  y  <  a,  b  <  y  <  00, 

4  & 

Vxy(0,y)  =y<7xy(0,y),  2/)  = /(?/)>  a  <  y  <  b, 

ayy(x,  0)  =  0,  crXy(xi  o)  =  0,  0  <  x  <  d, 

f  crxx(0,y)dy=  -  P, 

J  a 

€yy(xi  iOO)  60j 


(4a, b) 
(5a, b) 

(6a, b) 

(7) 

(8) 
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where  the  known  function  f(y)  defines  the  stamp  profile.  Note  that,  in  addition  to  f{y) 
the  external  loads  are  described  by  the  resultant  force  P,  the  remote  strain  e0  and  the 
crack  surface  tractions  given  by  (6a, b).  We  also  observe  that  the  unknown  functions  of 


the  problem  may  be  identified  as  follows, 

k  +  ldx(^X,°+)  V(X’0  ))  = 

0  <  x  <  d, 

(9a) 

k  +  1  <9rc(U(X,0+)  ))  =  ^(x)’ 

0  <  x  <  d, 

(9b) 

^xx(0,  y)  =  fz(y),  a  <  y  <b. 

(9c) 

In  the  following  sections,  we  will  derive  the  expressions  for  the  stresses  and 
displacements  in  the  terms  of  the  unknown  functions  fj,  (j  =1,2,3).  The  sum  of  the 
expressions  obtained  for  each  fj  must  satisfy  the  boundary  conditions  of  the  problem 
given  by  (4)-(8). 


2.1  The  contact  problem  (/i  =  0,  /2  =  0) 

In  this  section,  we  will  determine  the  stress  and  displacement  field  due  to  stamp 
loading,  namely  fz(y).  This  can  be  accomplished  by  using  Fourier  transforms.  The 
displacement  components  can  be  expressed  as 

1  f°° 

uz(x,y)  =  —  Uz(x,  p)exp(ipy)dp,  (10a) 

Tj-oo 

1  f°° 

vz{x,y)  =  —  /  Vz(x,p)exp(ipy)dp.  (10b) 

2  TTJ-co 

In  (10)  subscript  3  stands  for  the  displacements  due  to  stamp  loading.  Substituting  (10)  in 
(3)  following  ordinary  differential  equations  are  obtained: 
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(k  +  1)^  +  t(k  +  l)^f  -  p2(«  -  1)0»  +  2>P^  +  7V(3  -  k)Vz  =  0,  (11a) 

2ip- ^  +  7V(k  -  l)tfs  +  (k  -  1)^  +  7(k  -  1)^  -  P2(k  +  l)Va  =  0.  (1  lb) 

Assuming  a  solution  of  the  form  exp(srr)  following  characteristic  equation  is  obtained, 


(s2  +7 s-  p2  -  i\p\S3)  (s2  +  7 s-  p2  +  i\p\8z )  =0,  S3  =  \  j\ 


'3- K 

K  +  1  ’ 


(12a,  b) 


Roots  of  the  characteristic  equation  are  given  by 


si  =  ~\l~  tjvV  +  4p2  +  4z|p|<53,  K(si)  <  0, 


S2  =  ~  ^7  -  ^  VV  +  4P2  -  4i|p|<$3,  ^(s2)  <  0, 


(13b) 


S3  =  -  ^7  +  ^  \/72  +  V  +  4i|p|<$3,  »(*)  >  0, 


S4  =  -  ^7  +  ^  \/72  +  4p2  -  4i|p|<53,  ^(s4)  >  0. 


(13d) 


The  displacement  components  u%  and  v3  then  can  be  written  as 


t,  V)  =  7T  /  y'Miexp(sia:  + 

2^7-00  j=i 


^  /*oo  2 

u3(x,  y)  =  Y^MjNjQxp^SjX  +  ipy)dp, 


(14b) 


where 


Nj(p) 


z((k  +  l)s2  +  7(k  +  l)sj  +  p2(l  -  ac)) 

p(2sj  +  7(3  -  k)) 


,  (j=  1,2, 3, 4). 


Using  (2),  stresses  and  displacement  derivative  can  be  expressed  as  follows: 


<?xx3  (*,V) 


+  ipNj(3  -  K,))Mjexp(sjX  +  ipy)dp,  (16a) 


^(s,  y)  =  -^7 [  >>(3  “  «)  +  *pWj(/c  +  l))Miexp(sjx  +  ipy)dp,  (16b) 

K  -  1  27T  J-oo^ 

Z*00 

cr^ysCa:,  y)  =  /  ]T(«P  +  iVjsj)MJexp(sjx  +  ipy)dp,  (16c) 

^  J  ~ OO 


^  2  A00  ^  ^ 

«-w3(a,  y)  =  7-  /  ipJ^M;exp(sjx  +  ipy)dp. 


(16d) 


Using  the  boundary  conditions  given  by  (4)  and  (5a),  we  can  write, 

/OO  2 

yfafr  +  1)  +  ipNj{ 3  -  K,))Mjexp(ipy)dp  = 

■OO  j=l 

=  (  a  <  y  <  b  (17a) 

[0,  -oo<y<a:b<y<oc 

tojYttP  +  N^M^ipy^p  =  | ^Mv)’  _ b<y<00  <17b> 


We  express  Mj(p),  (j  =  1, 2)  in  the  following  form, 
Mj(p)  =  —  <M/>)  [  /3(t)exp(  -  ipt)dt. 

Po  Ja 

Then,  ipj(p),  (j  =  1, 2)  is  determined  from 

2 

+  1)  +  ipiVi(3  -  «))V»i(p)  =  («  ~  1). 
j=l 

2 

J^(ip  +  Njsj)ipj(p)  =  rj. 

3=  1 


(19b) 


Stresses  and  displacement  derivative  for  the  stamp  loading  can  now  be  obtained  using 


2.2  The  opening  mode  problem  (/2  =  0,  /3  =  0) 


In  this  section  we  will  determine  the  stresses  and  displacement  derivatives  due  to 
relative  displacement  derivative  of  the  crack  faces  in  y  -  direction,  namely  fi(x).  First, 
we  will  derive  the  expressions  for  stresses  and  displacement  derivative  for  a  crack  in 
infinite  plane.  Then,  the  solution  for  the  half-plane  ( x  >  0),  will  be  superimposed  to 
satisfy  the  boundary  conditions  at  the  free  surface  x  =  0.  In  the  solution  of  the  half-plane 
problem,  we  will  also  consider  the  symmetry  about  a>axis.  For  a  crack  in  infinite  plane, 
displacement  components  can  be  expressed  using  Fourier  transformations  as  follows: 

i  r°° 

u^\x,y)  =  —  Uil\uj,y)exp(iux)du,  (20a) 

27ry_oo 


Or,  y)  =  —  y  vit]  (w>  2/)ex p{iux)du, 


(20b) 


where  subscript  1  stands  for  the  opening  mode  problem  (i.e.,  fi(x)  ^  0,  f2{x)  =  0, 
fe(x)  —  0)  and  superscript  (i)  stands  for  the  infinite  plane  problem.  Substituting  (20)  in 


(11)  following  differential  equations  are  obtained: 

(k  —  1)  ^2  "b  (K  +  1)(t ~  +  (2 iu  +  7(3  —  ft))- —  =  0, 


(21a) 


(2 iu  +  7(ft  -  1))- 


(ft  +  1  )~tt  +  (K  ~  -  u?)Vi  =  0. 


(21b) 


Assuming  a  solution  of  the  form  exp  (ny)  the  characteristic  equation  is  found  to  be 
(n2  —  8\n  +  iu{^  +  iu))  ( n 2  +  Sin  +  iu(j  +  iu))  =  0,  Si  =  ^  •  (22a, b) 


Roots  of  the  characteristic  equation  are  given  by 


711  =  -  +  \\fkJ  -  4 iu^+Sj,  3? (ni)  >  0, 
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n2  =  +  \J 4w2  —  ~4iur/  +  6% ,  5R(n2)  >  0,  (23b) 

m  —  -  ^  J 4c ~2  -  4io;7  +  8\,  ^(^3)  <  0,  (23c) 

77,4  —  —  —  \  / 4u^  —  Aiury  +  <5^ ,  3^(714)  <  0.  (23d) 

^  Z  ’ 

Then,  for  y  <  0  and  y  >  0  stresses  and  displacements  can  be  expressed  as  follows: 


y  <0 

«P(*,  y)  =  —  ^Cjexp(nir/  +  iwx)du>, 

Z7T  J _ 00  j—i 

(24a) 

1  z*00  ^  A 

t>P (*»  y)  =  7T  YlC3AiQXV(niy  +  lux)dbJ-> 

^  J  —  OO  j— 1 

(24b) 

<&i(*>y)  =  ^  2^f  ^ytU^Mnjy  +  iux)du, 

(24c) 

°yyi(xi  y )  =  ^  2nJ  JZ,Syyi^QX^y  +  zu;x)du;’ 

(24d) 

1  f°° 

417  (*,  y)  =  M*)  7T  +  ^j)C,exp(niy  +  iu;x)da;, 

^  J  -OO  j— 1 

(24e) 

9  1  Z*00  A 

— u[r }  (x,  y)  =  —J  ^Ciniexp(njy  +  io?x)du;, 

(24f) 

2 

4xij  ~  yi(^(K  +  1)  +  A3rij(3  —  K,))Cj, 

3= 1 

(24g) 

4/ij  =  ]P(i^(3  -  «)  +  +  l))Cj, 

(24h) 

j=i 


where  superscript  (i  )  stands  for  y  <  0. 
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+  2  f°°  4  A 

u[l+)(x,  y)  =  —  /  '^2cjexp(njy  +  iux)du, 

Z7rJ-oo  j- 3 

(25a) 

+  ^  r°° 4  A 

V{f  \x,  y)  =  —  /  ^CJA;exp(niy  + 

^  J -co  j- 3 

(25b) 

y )  =  ^ 2irJ  S5Li;Mexp(^y  +  ««*0dw> 

(25c) 

4/A*’  v)  =  *  (_  1  27r/ooE'S^il(a;)exP(n^  + 

(25  d) 

^iv?  (*.  y)  =  /*(*)  X-  /  +  ^)QexP(^y  +  iux)du;, 

J-ooj =3 

(25  e) 

^-uf+)(x,  y)  =  -^~  f  'S2cjnjexp(njy  +  iu;x)dcu, 
dy  2t tJ-oo-^3 

(25f) 

=  y^(za)(K  +  1)  +  -  K))Cj, 

i=3 

(25g) 

^yii  =  X/^za;^3  -  «)  +  AjUj{K  +  l))Cj, 
j=3 

(25h) 

where  superscript  (i+)  stands  for  y  >  0.  In  (24)  and  (25)  Cj  ( j  = 

constants  and  A,  is  given  by 

1, 2, 3, 4)  are  unknown 

,  ,  _  n^(K- l)  +  (iu;7-u;2)(K  +  l) 

J  U  rij(2iu  +  7(3  —  k)) 

(26) 

For  a  crack  in  infinite  plane  and  for  the  opening  mode  following  boundary  conditions 
must  be  satisfied: 


—  00  <  X  <  oo, 


(27b) 


ax2(x,  °)  =  °) 


7Tit(u'')(x’0)-u'~){x'0))  =  0’ 


2/^0  d  /  (i+) 
K  +  1  02  V 


;(x,  0)  —  v[  \x,0 


—  OO  <  X  <  oo 


oo  <  £  <  0,  d  <  x  <  oo 


fi(x),  0  <  x  <  d 


(27c) 


(27d) 


We  first  express  the  unknown  constants  Cj(u>)  (j  —  1, 2, 3, 4)  as 

Cj(u)  =  f  /i(^)exp(  -  iurf)di. 

Jo 


(28) 


Then,  by  using  (24),  (25)  and  (27)  the  following  equations  can  be  obtained  to  determine 

4  2 

^(iu;(3  —  k)  +  Ajnj(l  +  k))Pj(u)  —  ^(iu;(3  —  /c)  +  -Ay^y(l  +  K))Pj(u )  =  0, 

;=3  j=i 


+  iwAj)Pj(u)  -  ^2(nj  +  iuAj)Pj(u )  =  0, 

j=3  J=1 


■ji’jM  -  £  V>>)}  =  !. 

3=  1  J 

io; jp4(u;)  +  Pz(u)  -  Pzfa)  -  P\(u) }  =  0. 


(29a) 

(29b) 

(29c) 

(29d) 


Using  (34),  (25)  and  (29)  the  stresses  and  the  displacement  derivative  for  the  opening 
mode  can  be  obtained.  Note  that,  if  we  solve  (29)  for  Pj(u)  then  substitute  in  (28)  and 
(24e),  (25e),  we  find  that 

<? («,  0)  =  &xyl  (*I  0)  =  0  (30) 
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which  is  expected  due  to  the  symmetry  about  x-axis.  For  the  opening  mode  problem 
tii(x,  y)  is  an  even  and  17  (x,  y)  is  an  odd  function  of  y.  In  order  to  satisfy  the  free  surface 
boundary  conditions  at  x  =  0,  we  will  now  superimpose  the  solution  for  the  half-plane 
(x  >  0)  on  the  solution  for  infinite  plane,  and  because  of  symmetry  in  the  half-plane 
(x  >  0)  solution  we  will  only  consider  y  >  0.  Hence,  displacement  components  for  the 
half-plane  can  be  written  by  using  the  following  Fourier  cosine  and  sine  integrals: 

poo 

u[h\x,y)  =  /  ulh\x,a)cos(ay)da,  (31a) 

Jo 

poo 

v^\x,  y)  =■  /  (x,  a)sin(ay)da,  (31b) 

Jo 


where  superscript  (h)  stands  for  the  half-plane  problem.  Substituting  (31)  in  (3)  following 
ordinary  differential  equations  are  obtained: 


,  ,.*u}k) 

(K  +  l)-lxr 


+  t(k 


-  a2  {k-  l)u[h]  +  2a-J—  +  7( 3  -  K)Vfh)  =  0, 


d2V{h) 


2a^  _  7tt(K  _  i)^  +  (K  .  +  7(K  _  i)-£. 


a2(K  +  l)Vi(/l)  =  0. 


(32b) 


Assuming  a  solution  of  the  form  exp(px)  we  obtain  the  characteristic  equation  as, 

( p 2  +  7 p  —  a2  —  (p2  +7 p  —  a2  +  ia6\ )  =  0,  (33) 


where  is  given  by  (22b).  Roots  of  the  characteristic  equation  are  found  to  be 


&(pi)  >  0, 

&(P2)  >  0, 


(34a) 

(34b) 
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Pz  =  -  ]p- \\fl2  +  4a2  +  4ia5i,  3fJ(ps)  <  0, 

P4  =  —  -7  —  -  \/72  +  4a2  —  4ia6i,  5ft(p4)  <  0. 

2  2 


(34d) 


The  stresses  and  displacements  for  the  half-plane  problem  are  then  expressed  as  follows: 

I 

poo 

uf\x,y)=  l  (B3exp(p3x)  +  54exp(p4z))cos(m/)da:,  (35a) 

Jo 

poo 

v(fl\x,y)  =  (BzD3exp(pzx)  +  B4DAexp(p4x))sm(ay)da,  (35b) 

Jo 

y)  =  -^T  [  +  1)Pi  +  Dia<<3  -  K))Bjcxp(pjx)cos{ay)da,  (35c) 

K-lJo  ^ 

<Jyy\ (x,  y)  =  [  B(3-  #0Pi  +  +  l))5iexp(pjx)cos(ay)dQ:,  (35d) 

^  ~  -L  JO  1=s 


=  /*(*) 


rco_ji 

/  -  a)-Bjexp(pj 

40  j=3 


?x)sin(ay)do:, 


poo 

4ft)(x,y)=  -  /  a(jB3exp(p3a:)  +  54exp(p4x))sin(ay)d£ 
Jo 


where  Bj  ()  =  3, 4)  are  unknown  constants  and  Dj  is  given  by, 


Dj  = 


P2j(k  +  1)  +  a2(  1  —  k)  +  7Pj(l  +  K) 
«(2pj  +  7(3  -  k)) 


For  x  >  0  and  y  >  0,  the  total  stress  and  displacement  fields  can  be  obtained  by  adding 
the  equations  (25)  and  (35),  that  is 


ui(x,y)  =  u(f\x,y)+uf\x,y), 


(37b) 


v\ (x,y)  =  vf\x,y)  +  v[h)(x,y), 
&xxi  (x,  y)  =  42  (x,  y)  +  critic  !/)» 
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(37c) 


(Jxyl  (x,  y)  =  cr^yi  (X,  y)  +  4i  (*>  l/)>  (37d) 

(*,  y)  =  42  +  4  (x’  y )  •  (37e) 

The  constants  -Bj(a),  (j  =  3, 4)  are  determined  by  using  the  free  surface  boundary 
conditions  as  follows, 

<r*xi(0, 2/)  =  42(0,  y)  +  4i(°>  y)  =  °>  0  <  y  <  oo,  (38a) 

Oxyl (0, y)  =  42 (°>  ^)  +  41  (0, y)  =  o,  0  <  y  <  00.  (38b) 


Note  that  due  to  symmetry  we  only  consider  0  <  y  <  oo.  Using  (25c),  (25e),  (35c),  and 
(35e)  and  after  simplifications  using  by  MAPLE,  (38)  is  reduced  to  following  form: 

X2(k  +  l)pi  +  £)iQ!(3  -  K))BAa)  = 

i=3 

1  K_i  W  r00 

= - ~  — — /  h(t)dt  Fxxi(u, <*)exp(  -  iut)du,  (39a) 

7T2  2/X0  Jo  J-oo 

4  I  I  /*<*  /■o° 

V (DjPj  -  a)Bj(a)  =  -  —  —  /  /i(t)dt  /  Fwi(w,  a)exp(  -  iut)du  (39b) 

^3  TT  ^/A)  Jo  J-oo 

where, 


Fxxii^,  a) 


AioPto 
D(uj,  a)  ’ 


-fxyl  (k*j  tu) 


4a>o:(u;  —  *7) 
•D(u>,  a) 


D(a;,  a)  =  uA  —  2?7u;3  +  (2o:2  —  72)^2  -  2ia2'fu  +  a4  +  a272 


3  —  n 
K  T  1 


(40a, b) 
(40c) 


The  inner  integrals  in  (39)  are  evaluated  in  closed  form  using  the  theory  of  residues  and 
(39)  is  reduced  to  following  form  to  determine  the  unknown  constants  Bj(a),  ( j  —  3, 4): 

4 

y2((«  +  1  )Pj  +  DjOt{ 3  -  n))B*(a,t )  =  Rxx\(a,t),  (41a) 

3= 3 
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(41b) 


V  'jjDjPj  a)Bj  (^s  “0  —  Rxyl(Oi,t). 

3= 3 

Bj(a)  is  now  defined  as 

Bj(a)  =  B*(a,t)ex] :p(Q  - 


(42a) 


-Rzxl  (^J  t) 


1/^  —  1  or 


K  +  1  AiA2(A^  +  A|) 


|7A2COs(A2i)  +  (7A1  —  2(A^  +  A2)  sin 


(42b) 


■E^xyl  (^?  t) 


2  1 


7T  K  +  1  AiA2(A^  +  A|) 

X  |a2(Ai  +  A2  +  72/4)cos(A2t)  -  Ai  (A2  +  \\  -  72/4)sin 


where 


R\  +  R2 


IR1-R2 


(43a, b) 


Ri  =  \j  (72/4  +  a2)2  +  ck272(3  —  k)/(ac  +  1), 


(43c) 


i?2  =  72/4  +  a2. 


(43d) 


This  completes  the  formulation  for  the  opening  mode  problem.  Stresses  and 
displacements  are  given  by  (37)  and  unknown  constants  are  given  by  (29)  and  (41). 

2.3  The  sliding  mode  problem  (/i  =  0,  /3  =  0) 


In  this  section,  we  will  determine  the  displacements  and  stresses  due  to  relative 
displacement  derivative  of  the  crack  faces  in  ^-direction,  namely  ^(x).  First  we  will 
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derive  the  expressions  for  stresses  and  displacements  for  a  crack  in  infinite  plane.  Then, 
the  solution  for  the  half-plane  (x  >  0)  will  be  superimposed  to  satisfy  the  boundary 
conditions  at  the  free  surface  x  =  0.  Again,  symmetry  will  be  considered  in  the  solution 
of  the  half-plane  ( x  >  0)  problem.  Following  a  similar  procedure  as  given  in  Section  2.2, 
stresses  and  displacements  for  the  infinite  plane  can  be  written  as  follows, 


y  <  0 


2  2  ^ 

uf\x,  y)  =  —  JZ^exp  (njy  +  iux)du, 

^  J  — 00  j— i 

(44a) 

v(2n(x,  y)  =  —  /  ^EjAfMnjy  +  iux)du, 

J  —  OO  j=l 

(44b) 

<r£d(x>  y )  =  1  2nJ  J2SU(^Mnjy  +  iwx)du, 

(44c) 

vm  ^y)  =  + 

(44d) 

2  f°°  2  a 

ofj (X,  y)  =  fi(x)  7T  /  +  ^Aj)EjQMnjy  +  iux)du, 

J  —  oo  j—i 

(44e) 

q  ^  roo  2  ^ 

— uf 5  (ar,  y)  =  —  /  V^nJexp(njy  +  ia;x)du;, 

9y  2tt  J  _qo 

(44f) 

2 

5^  =  +  1)  +  Ajrij(3  -  n))Ej, 

3= 1 

(44g) 

”  ^^(^(3  “  k)  +  i4j7lj(/C  + 

(44h) 

j=i 

where  subscript  2  stands  for  the  sliding  mode  or  mode  II  problem  and  superscript  (i~) 
stands  for  y  <  0. 
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y  >  o 


+  ^  roc  4  ^ 

uf\x,  y)  =  —  ^EjexpinjV  +  ivx)du, 

^  J  -co  j— 3 

(45a) 

V(p(x,  y)  =  —  J2EJAJeMnjy  +  iux)du, 

J- 00  j=3 

(45b) 

<4*2  (*,  y)  =  J  2„J  J2SU(“)QMnjy  +  *w*)dw, 

(45c) 

4/2  (x’  y)  =  \  27r/oo54544(a,)exP(n^  +  ^)du;, 

(45d) 

+  2  r00  4  ^ 

&xy2  (*»  y )  =  M*) 7T  /  54 (nJ  +  iuAj)Ejex]j(njy  +  iwx)du;, 

^J-ooj=3 

(45e) 

5  T  1  r00  ^ 

—u[p{x,  y)  =  —  V£inJexp(n;y  +  iwx)dw, 

2tt  j_0O  ^._3 

(45  f) 

Exx2j  =  54(Mk  +  1)  +  AjTlj{Z  -  k))-E), 
i=3  • 

(45g) 

4 

Eyy2j  ~  V  ](^(3  —  k)  +  AjTlj^K  +  l))£^’, 

(45h) 

j=3 


where  superscript  (i+)  stands  for  the  half-plane  (y  >  0),  rij,  (j  =  1, 2, 3, 4)  is  given  by 
(23)  and  Aj  is  given  by  (26).  The  unknown  constants  Ej,  ( j  =  1,  2, 3, 4)  axe  determined 
using  the  following  boundary  conditions: 


<4£?(*.  °) = 44  (*,  o), 

—  OO  <  X  <  00, 

(46a) 

=a,Sf2(**°)i 

—  00  <  X  <  00, 

(46b) 

>(x,0)  ^'(l’0); 

II 

—  oo<x<0,  d  <  x  <  oo 

0  <  x  <  d 

(46c) 
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which  is  expected  due  to  the  symmetry  about  the  x-axis.  For  the  sliding  mode  problem 

u2(x,  y)  is  an  odd  and  v2(x,  y)  is  an  even  function  of  y.  In  order  to  satisfy  the  free  surface 

boundary  conditions  at  x  =  0,  we  will  now  superimpose  the  solution  for  the  half-plane 

( x  >  0)  on  the  solution  for  the  infinite  plane,  and  because  of  symmetry,  in  the  half-plane 

(x  >  0)  solution  we  will  only  consider  y  >  0.  Hence,  displacement  components  for  the 

half-plane  can  be  written  using  Fourier  sine  and  cosine  integrals  as  follows: 
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where  superscript  (h)  stands  for  the  half-plane  problem.  Substituting  (50)  in  (3)  following 
ordinary  differential  equations  are  obtained 


(«  +  l)^f  +  +  -  D^»  -  2^  -  7q(3  -  >  -  0, 


2^  +  7Q(K  _  1)vm  +  (K  _  ^  +7(K  _  d-SlI  -  «»(«  +  1}vf  >  -  o. 


(51b) 


Assuming  a  solution  of  the  form  exp(ta),  the  characteristic  equation  is  obtained  as 


(t2  +  jt  —  o?  —  za^i)  ( t 2  +7 1  —  a2  +  ia6\ )  =  0, 


where  6i  is  given  by  (22b).  Roots  of  the  characteristic  equation  are  obtained  as  follows 


ti  =  —  ^7+  \  yfl 2  +  4a2  +  Aia6\ ,  3ft(ti)  >  0, 

z  z 

t2  =  —  x7+  ;;  VV  +  4a2  -  4za<5i ,  Kfo)  >  0, 

z  z 

*3  =  -  ^7  -  7;  vV  +  4o;2  +  4*0^1,  R(t3)  <  0, 

t4  =  —  ^7  —  ^  \/72  +  4a2  —  4ia6i ,  3£(f4)  <  0. 

z  z 


(53b) 


(53d) 


For  the  half-plane  solution  the  stresses  and  displacements  can  now  be  expressed  in  the 
following  form 


roo 

\x,y)=  /  (G3exp(t3x)  +  G4exTp(Ux))sm(ay)da, 
Jo 


poo 

vf\x,y)  =  /  (G3#3exp(f3x)  +  G4tf4exp(t4x))da:, 

Jo 

<4*2^,  y)  =  f  Y\{{K  +  l)*j  -  <*#,(3  -  i^))Gjexp(tjx)sin(ay)da, 
K~  LJo  ,=3 


j(x,  y)  =  £3^  ]T((3  “  K)ti  ”  a#i(K  +  1 ) ) ^ jexP (^j^) sin( ay) da, 


<4y2(*>  2/)  —  ^(x  )  E(“+^<  j)Gjexp(^x)cos(ay)do;, 

•/O  j=3 

5  A00 

—u£\x,y)  =  /  a(G3exp(t3x)  +  G4exp(f40)cos(m/)da, 


where  G7,  (j  =  3, 4)  are  unknown  constants  and  Hj  is  given  by, 

rr  ,  x  7*i(«  +  1)  +  a2(l  -  *0  +  *?(«  +  !)  /  -o.N 

—  7—  ■  rr  rr  ?  U  ~ 

J  a(2ij  +  7(3  -  «)) 


(54b) 

(54c) 

(54d) 

(54e) 

(54f) 

(55) 


For  x  >  0  and  y  >  0,  the  total  stress  and  displacement  fields  can  now  be  obtained  by 
adding  (45)  and  (54)  as 


u2(x,  y)  =  uf\x,  y)  +  u}\x,  y),  (56a) 

v2(x,  y)  =  vf\x,  y)  4-  vf\x,  y),  (56b) 

°xx2{x,  y)  =  a  fxl (x,  y)  +  a{x hJ2(x,  y),  (56c) 

oxy2  (x,y)  =  <JXy2  (x,y)  +  o%  (x,y),  (56d) 

<Xyy2(x,  y)  =  cr^  (x,  y)  +  cr^(x,  y).  (56e) 


The  constants  Gj(a),  (j  —  3,4)  are  determined  by  using  the  free  surface  boundary 
conditions  as  follows: 
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<7**2  (o,  y )  =  42  (°>  y)  +  4*2  (°» y )  =  °> 


0  <  y  <  oo, 


(57a) 


<7^2  (o,  ?/)  =  42  (°>  y) +  42  (°» y)  =  °>  0  <  v  <  °°- 


(57b) 


Note  that  due  to  symmetry  we  only  consider  0  <  y  <  oo.  Using  (45c),  (45e),  (54c)  and 
(54e)  and  after  simplifications  by  using  MAPLE,  (57)  is  reduced  to  following  form: 

£((*  +  1)*;  -  Hja{Z  -  «))<?»  = 

i= 3 

=  f2(t)dt  [  Fxx2(uj,a)exp(-iujt)du,  (58a) 

2/i0  Jo  7-oo 

4  11/*^  r°° 

+  a)Gj-(a)  =  -  ^  —  J  hit)dt J  Fxy2{u,a)exv{-iut)du,  (58b) 


where, 


Fxx2(u,a)  = 


D(u>,  a)  ’ 


Fxy2(ijJ,  a)  = 


Aa2(iu  +  7) 
D(w,  a)  ’ 


(59a, b) 


and  D(u,  a)  is  given  by  (40c).  The  inner  integrals  in  (58)  are  evaluated  in  closed  form  by 
using  the  theory  of  residues  and  (58)  is  reduced  to  following  form  to  determine  the 
unknown  constants  Gj(a),  (j  =  3, 4): 

£((k  +  l)tj  -  Hja( 3  -  K))G*(a,t)  =  Rxx2(oc,  t),  (60a) 

j= 3 


+  a)G*(a,  t )  =  Rxy2(a,  t). 
j= 3 


(60b) 


where 


G»  =  ^  f  G|(a,()ex p((7  -  A .i)i)/j(t)<ft 


(61a) 


R„2(«,t)  =  -  +  +  AismM}._  (61b) 

„  ,  .A  11  O? 

Rxy2{a,  t)  -  ^  k  +  x  Ai  As(A2  +  A2)  X 

x  |7A2cos(A2t)  +  (2(Af  +  Xf)  +  7A1  )sin(A2i) J.  (61c) 

In  (61)  Ai  and  A2  are  given  by  (43).  This  completes  the  formulation  for  the  sliding  mode 
problem.  Stresses  and  displacements  are  given  by  (56)  and  unknown  constants  are  given 
by  (48)  and  (60). 

3.  Derivation  of  the  singular  integral  equations 


The  stress  and  displacement  fields  for  the  contact,  opening  mode  and  sliding  mode 
crack  problems  are  given  in  Section  2  in  terms  of  the  unknown  functions  f2{t)  and 
fz(t).  The  total  stress  and  displacement  fields  can  now  be  expressed  as  follows. 


u(x,  y)  =  ui(x,  y)  +  u2(x,  y)  +  u3(x,  y), 

0  <  x  <  oo, 

0  <  y  <  oo, 

(62a) 

v(x,y)  =  vi(x,y)  +  v2(x,y)  +  v3{x,y), 

0  <  x  <  oo, 

0  <  y  <  oo, 

(62b) 

&xx{xi  y )  =  <Xxxi(X) y )  4"  (Xxx2.(.xi  y)  ~b  Gxxafai y)i 

0  <  x  <  oo, 

0  <  y  <  oo, 

(62c) 

&xy(x>y)  =  Gxyl(Xi  y)  +  crxy2(^>J/)  "I"  O'xy^X,  y), 

0  <  x  <  oo, 

0  <  y  <  oo, 

(62d) 

&yy(xi  y)  —  <7j/j/l(a;5  y)  ~b  <xyy2{X)  y)  "b  GyyziXi  2/)  "b  JE(x)Ci o, 

0  <  x  <  oo, 

0  <  y  <  oo, 

(62e) 

where  expressions  due  to  fi,  /2  and  /3  are  given  by  (37),  (56)  and  (16)  respectively.  In 
this  section,  the  problem  will  be  reduced  to  three  singular  integral  equations  using  the 
boundary  conditions  (5b)  and  (6),  i.e.,  by  using. 


<7yy(x,0)  =  0,  0  <x<d,  (63a) 

axy(x,  0)  =  0,  0  <  x  <  d,  (63b) 

■^rrS-«(0,  y)  =  f(y),  a<y<b.  <63c) 

k  +  1  ay 

Considering  the  expressions  given  in  Section  2,  stresses  and  displacement  derivative  can 
be  written  in  the  following  general  form, 

<jyy{x,y)  -  f  j(x,y,t)fj(t)dt+  f  hz{x,y,t)fz{i)dt  +  E(x)e0,  (64a) 

J 0  j=\  J  & 

axy(x,y)  =  [  y,  t)fj(t)dt  +  f  k2z{x,y,t)fz(t)dt,  (64b) 

JO  j—\  J & 

-^-S-u(x,y)=  f  y2kzj(x,y,t)fj(t)dt+  f  kzz(x,y,t)fz(t)dt.  (64c) 

1  +  ndy  Jo  pi  Ja 

Then,  using  (63)  integral  equations  can  be  expressed  as  follows: 

ayy(x,  0)  =  f  'Vklj(x,0,t)fj(t)dt+  f  k13{x,0,t)fz(t)dt  +  E(x)e0  =  0, 

J 0  Ja 

0  <  x  <  d,  (65a) 

rd  2  rb 

axy(x,0)  =  /  YM*>  t)fj(t)dt  +  /  k2z(x,0,t)fz{t)dt  =  0, 

Jo  Ja 

0  <  x  <  d,  (65b) 

4fi2-^-«(0,v)=  fdf'hj(o,y,t)mdt+  /  MO,  »,*)/»(*)<«  =  /(»). 
l  +  yo  •/« 

a  <  y  <  b.  (65c) 
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In  the  previous  section  it  was  shown  that  because  of  symmetry  oyy2 (x ,  0)  =  0,  for 
fi(t)  =  0,  f2(t)  ±  0,  /3(f)  =  0  and  a*yi(x,0)  =  0  for  /i(i)  ±  0,  f2(t)  =  0,  fz(t)  =  0, 
from  which  it  follows  that 


kn(x,0,t)  =  0,  k2l(x,0,t)  =  0.  (66a, b) 

The  expressions  for  the  other  kernels  kij(x,  y,  t)  (kz](x.  0,  t )  for  i  =  1,2,  and  fcjj( 0,  y,  t) 
for  i  =  3)  will  be  given  in  the  following  sections.  The  kernels  will  be  expressed  in  the 
following  general  form, 

rco 

fcii(x,0,t)=/  Kij(x,  t,  p)dp,  (i  =  1,2,  i  =  1,2,3),  (67a) 

JO 

roo 

kij( 0,  y,t)=  /  Kij(y,  t,  p)dp,  (t  =  3,  j  =  1, 2, 3).  (67b) 

With  the  exception  of  one  case  (that  being  £133(0,  y,t)  for  7  <  0)  the  integrands  in  (67) 
are  bounded  and  continuous  for  p  <  00  and  integrable  at  p  =  0.  The  singular  nature  of 
the  kernels  ktj  is  therefore  determined  by  examining  the  asymptotic  behavior  of  the 
integrands  as  p  tends  to  infinity.  In  the  following  sections  we  will  give  the  expressions 
and  details  of  the  asymptotic  analyses  for  the  kernels. 

3.1  ku(x,y,t) 

We  first  express  kn(x,  y,  t)  as  follows 

kn(x,y,t)  =  k$(x,y,t)  +  k[ \\x,y,t),  (68) 

where  k^  is  obtained  from  the  infinite  plane  solution  and  k[ ^  is  obtained  from  the  half¬ 
plane  (x  >  0)  solution.  Referring  to  (25d),  fcJJ  (x,  y,  t)  can  be  written  as 
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K{{x,y,t) 


_  k  +  1  exp  (7a:)  /'°°  « 


K  —  1  Ak 


/°° 

y)Q\p(iu(x  -  t))du, 

•00 


<f>n  (u>,  y)  =  3  -  k)  +  AjUj(  1  +  At))P;(u;)exp(niy),  (69b) 

7=3 

where  n3,  ( j  -  1,2, 3, 4),  A3  and  P3  are  given  by  (23),  (26)  and  (29),  respectively. 
Changing  the  limits  of  integration  (69a)  can  be  written  as, 


k[ii(x,y,t) 


k  +  1  exp  (7a;) 
k  —  1  4n 


j  Pn1(w,y)cos(o;(a:  -  t ))  + 

+  K$2(u,  y)sin(u(x  -  t))^du, 


where. 


Kni(u,  y)  =  4>n(u,  y)  +  </>§(  -  w,  y), 

y )  =  *  (^11  (w»  y)  -  tfi  ( -  y))  • 


(71b) 


In  order  to  extract  the  singular  terms  we  expand  and  into  series  as  to  —>  co. 
Following  asymptotic  expansions  are  obtained  by  using  MAPLE: 

!/)=(-  +  ^  +  -+  ^+  °(A)  W  -  “’ll).  <72a) 


cj  or 


tfSfW)  =  /20  +  if  +  •••  +  +  O 


exp(  —  coy), 


where  the  leading  terms  are 


/20  = 


4(k~1) 

K  +  1  1 


2(k  —  1)7 
K  +  1 


(73a, b) 


Subtracting  the  asymptotic  expansions  from  the  integrands  in  (70),  using  integration 
cutoff  points  for  the  infinite  integrals,  evaluating  some  of  the  integrals  in  closed  form  and 
taking  the  limit  as  y  — ►  0,  after  some  manipulations  (x,  y,  t )  is  reduced  to 


where  A^2  is  another  integration  cutoff  point,  70  is  the  Euler  number,  [10].  Second 
integrals  in  (75  a)  and  (75b)  will  be  neglected  in  numerical  computation  for  sufficiently 
large  values  of  integration  cutoff  points.  Third  integrals  in  (75a)  and  (75b)  are  evaluated 
in  closed  form.  The  expressions  used  in  the  evaluation  of  these  integrals  are  given  in 
Appendix  A.  Referring  to  (35d)  and  (42a),  k[\\x,  y,  t)  can  be  written  as  follows 


k  +  1  exp(7x)  /■“  (A) 

K-l  2  io  11 


(a,t,x)cos(ay)da 


(76) 
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where 


K^\a,t,x)  =  (a,  t)exp(pjX  +  (7/2  -  Ax )t), 

j= 3 


4>w)(a)  =  Pifi  ~  K)  +  Dja(l  +  k), 


0  =  3,4). 


(77b) 


Pj,  (j  =  1,2, 3, 4),  Dj{ a)  and  B*{a,t)  are  given  by  (34),  (36)  and  (41)  respectively.  In 
order  to  extract  the  singular  terms  we  expand  x)  into  series  as  a  tends  to 

infinity  as  follows: 

tfn)o°(M,aO  =  |h2a2  +  ^a  +  /i0  +  ^-  +  ^  +  ---  +  ^7  +  °^)|  x 


x  exp(7(t  —  x)/2  —  (t  +  x)a), 


where  the  coefficients  of  the  expansion  are  also  fimctions  of  x  and  t.  The  leading  term  is 
given  by, 

_  16  k  -  1  sin(6ia:/2)sin(6xt/2) 

2  ~~  7k  +  1  8\ 

61  is  given  by  (22b).  Subtracting  K$)o°  from  the  integrand  in  (76),  evaluating  some  of 
the  integrals  in  closed  form,  using  an  integration  cutoff  point  for  the  infinite  integral 
and  taking  the  limit  as  y  — »  0,  after  some  manipulations  can  be  written  as 

*8>(«,o,*)  =  +  4?W)  - 


-  hiexp(y(t  -  i)/2)Ei(  -  +  xlj  |, 


where  Ei()  is  the  exponential  integral  [10]  and, 


Hi  («, t) 


(t  +  xf  ( t  +  X )2  t  +  : 


exp (7(*  ~  *)/2), 


(81a) 
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rAT 

J$\x,t)  =  K['l)*(cx,t,x)da 

Jo 

r°°  r  1 

+  /  ftfnHM*®)  -  #ff)o°(<Mj®) 

•M?  l  j 

POO 

+  /  r^°°(a,t,x)da,  (81b) 

./<> 

.£^*(0:,  t,  x)  =  i^{f}(Q;,t,x)  -  (/i^cc2  +  h\a  +  /i0)exp(7(i  —  x)/2  —  (£  +  x)a), 

(82a) 

/^“(a,  t,x)  =  l- |  +  ^|  +  •  •  •  4-  ^ )exp(7(t  -  x)/2  -  (t  +  x)a).  (82b) 

X1  ( cr  cr  a'  J 

Second  integral  in  (81b)  will  be  neglected  in  numerical  computation  for  a  sufficiently 
large  value  of  A^.  Third  integral  in  (81b)  is  evaluated  in  closed  form.  The  expression 
used  in  the  evaluation  of  this  integral  is  given  in  Appendix  A.  k^(x,  0,  t)  and  k[^  (x,  0,  t) 
are  given  by  equations  (74)  and  (80)  respectively.  Adding  these  two  equations  A:n(x,  0,  t) 
may  be  expressed  as, 

fcn(x,0,t)  =  exp(7x) |  —  — - h  hUs(x,t)  +  huf(x,t)  i,  (83) 

^  7T  t  —  X  J 

where, 

hiu(x,t)  =  ^^rn(x,t),  (84a) 

huf(x,  t)  =  —  —  j  —  /liln^inl^  —  +  Jn\(xA)  +  Jin(xi  j"  + 

+  {■/n)(x,t)  -  h!exp(7(t  -  x)/2)Ei(  -  A$(t  +  x))  (84b) 


28 


The  terms  in  (84)  are  given  explicitly  by  (75)  and  (81).  Note  that  the  first  term  in  (83)  is 
the  Cauchy  singularity  associated  with  a  crack  in  infinite  plane,  as  for  the  second  term  we 
can  write, 


lim  hiis(x,t) 

(x,t)—*0 


if  1  2 1 _ 4t2  \ 

iryt  +  x  +  (t  +  x)2  (t  +  xf) 


(85) 


This  term  becomes  singular  as  x  and  t  simultaneously  approach  zero,  and  is  the  standard 
expression  found  for  edge  cracks  in  homogeneous  materials,  (see  for  example,  equation 
(23a)  in  Dag  and  Erdogan  [8]). 


3.2  k13(x,y,  t ) 


Referring  to  (16b),  fc13(®,  y,  t)  can  be  written  as 
ki3(x,  y,  t)  =  eXP^~-  f  0i3 (p,  x)exp(ip(y  -  t))dp,  (86a) 

«  —  1  27 r  J-oo 

2 

013 (p,  x)  =  ^(ipNjin  +  1)  +  Sj( 3  -  K,))ipj(p)exip(sjx),  (86b) 

i= i 

where  s3  (j  =  1,2, 3, 4),  Nj  and  ip3{p)  are  given  by  (13),  (15)  and  (19),  respectively. 
Changing  the  limits  of  integration  in  (86a)  and  taking  the  limit  as  y  — >  0,  and  rearranging, 
kiz(x,  0,  t)  can  be  written  as 

k13(x,0,t)  =  ^ {^i3i(p,3;)cos(pt)-it:132(p,x)sin(pt)}dp,  (87) 

where, 

#131  (p,  X)  =  01 3(P,  X)  +  013 (  -  P,  *),  (88a) 

#132 (p> z)  =  i(0i3(P)  a:)  -  0i3(  -  P,  z))-  (88b) 
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In  order  to  extract  the  singular  terms  we  expand  Km(p, x)  and  K^ip, x)  into  series  as 
p  — >  oo.  Following  asymptotic  expansions  are  obtained  by  using  MAPLE, 


(93b) 


(93c) 


(93d) 

where  Ai32  is  another  integration  cutoff  point  and  the  remaining  terms  are  given  by, 


#1*31  (p,  x)  =  #131  (p,  x)  -  (dnP  +  rfio)exp(  -  (7/2  +  p)x),  (94a) 

r$M*)  =  {  7  +  7  +  •  •  •  +  757  }«p(  -  (t/2  +  »)*).  <94b> 

#132 (p>  *)  =  #i32(p,  x)  -  (^2iP  +  d20  +  <Wp)exp(  -  (7/2  +  p)x),  (94c) 

r&M  =  {  7  +  7  + ' ' '  +  77  }exP<  “  (t/2 + />)*)•  <94d> 


Second  integrals  in  (93c)  and  (93d)  will  be  neglected  in  numerical  computation  for 
sufficiently  large  values  of  A131  and  Ai32.  Third  integrals  in  (93c)  and  (93d)  are  evaluated 
in  closed  form.  The  expressions  used  in  the  evaluation  of  these  integrals  are  given  in 
Appendix  A.  Also  note  that. 
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lim  hns(x,  t)  =  -<  -  -  ,2~2  ~V-f  2  2Vi 

(x,t)-»0  7T  [(x2  +t2)  {X1  +  t2) 


This  term  becomes  singular  as  x  and  t  simultaneously  approach  zero,  and  is  the  standard 
expression  obtained  for  the  homogeneous  materials  (see,  for  example,  equation  (23b)  in 
Dag  and  Erdogan  [8]). 

3.3  k22(x,y,i ) 

We  first  express  k22  {x,  y,  t)  as  follows 


k22(x,  y,  t)  =  kg (x,  y,  t )  +  k\2\x,  y,  t),  (96) 

where  kyl  is  obtained  from  the  infinite  plane  solution  and  k^2  is  obtained  from  the  half¬ 


plane  (x  >  0)  solution.  Referring  to  (45e)  k%(x,  y,  t)  can  be  written  as 
kjl (x,  y,  t)  =  (k  +  1) £X^J-  f  <t>22(u’y)exV(ioj(x  ~  ’ 


(97a) 


\(oJ,y)  =  J2(nj  +  iwAj  )Zj(u)  exp  (rijy), 

j=3 


(97b) 


where  nj}  (j=  1,2, 3, 4),  A3  and  Zj  are  given  by  (23),  (26)  and  (48),  respectively. 
Changing  the  limits  of  integration  (97a)  can  be  written  as, 


kj§(x,y,t)  =  (k  +  i)exP|lfl^  {  K^(u,y) cos(u(x  -t))  + 

+  K$2(u,  y)s\n(u(x  -  t))|dw, 


where. 


i^22i  (w’  y) =  02  (w>  y)  +  (-u>y)’ 
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#222  y)  =  *  {$22  {u,y)~  $22  (  “  2/) )  •  (99b) 

In  order  to  extract  the  singular  terms  we  expand  K^i  and  into  series  as  u  — »  oo. 
Following  asymptotic  expansions  are  obtained  by  using  MAPLE 

("■ »)  =  { kf  +  k~j  +  •  •  +  £rf  +  0  (js )  }«p<  - «-»).  <100a> 

*t>. P)  =  {^o  +  ^  +  '  ■  ■  +  ^  +  o(^)  }exp(  -  (100b) 


where  the  leading  terms  are 


^20  =  - 


4 

«  +  1  ’ 


ku  — 


27 

«  +  1 


(101a, b) 


Subtracting  the  asymptotic  expansions  from  the  integrands  in  (98),  using  integration 
cutoff  points  for  the  infinite  integrals,  evaluating  some  integrals  in  closed  form  and  taking 
the  limit  as  y  — >  0,  and  after  some  manipulations  k^l  (z,  y,  t)  is  reduced  to, 


1  1 
7 Tt  —  X 


l  r 


r(0 


(0 


+  (k  +  1)—  |^22l(a:>^)  +  ^222^X1^)  ~ 


—  fcnln 


(102) 


where, 


J2(2)1(x,f)=  /  #22l(^°)cOsMx-t))dw 


4« 

221  (0 


0 


poo  r  -| 

+  #221  °)  -  #221°°(^  °)  cos(u;(x  -  t))dw 

JA{1  l  j 

#22i°°(u;)  °)  ~  kn/oj  cos (v(x  -  t))du 
f  fA™  cos  (a)  —  1 

-MW.  ~ 


+ 


poo 

JA® 

J  A221 


-da  k 
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(103a) 


A 22  OM)  =  f 

Jo 


!(0 

*222 


jF^22 (w>  °)  ~  ^20  sin(u;(x  -t))du 


roo  - 

+  /„  irg2(w,0)-ifS”(a,,0)  sin (u(x-t))d 
J4W  L  J 

«/  JI999 

+L 


222 

00 

(0  L 

222 


^22°° (wj  0)  -  ^20]  sin(o;(x  -  t))du, 


(103b) 


A221  and  A222  are  integration  cutoff  points  and  70  is  the  Euler  number,  [10].  Second 
integrals  in  (103a)  and  (103b)  will  be  neglected  in  numerical  computation  for  sufficiently 
large  values  of  integration  cutoff  points.  Third  integrals  in  (103  a)  and  (103b)  are 
evaluated  in  closed  form.  The  expressions  used  in  the  evaluation  of  these  integrals  are 
given  in  Appendix  A.  Referring  to  (54e)  and  (61a),  k^(x,  y,  t )  can  be  written  as  follows 

k^(x,y,t)  =  («+  1  )exp^7X?^  K^\a,t,x)cos(ay)da,  (104) 


where 

4 

K${a,  t, x)  =  5^(a  +  Hjtj)G*(a,  t)exp(i,x  +  (7/2  -  Ai)t),  (105a) 

j= 3 


tj,  (j  =  1, 2, 3, 4),  Hj(a )  and  G*(a,  t )  are  given  by  (53),  (55)  and  (60),  respectively.  In 
order  to  extract  the  singular  terms  we  expand  K^idjtjX)  into  series  as  a  tends  to 
infinity  as  follows: 


K^(c,t,x)  = 


*  9  *  ml 

m2or  +  mxa  +  mo  4-  — 


+  ^j  + 

or 


+ 


m.7 

If 


+  0 


x  exp(7(f  —  x)/2  —  (t  +  x)a),  (106) 


where  the  coefficients  of  the  expansion  are  also  functions  of  x  and  t.  The  leading  term  is 
given  by, 
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(107) 


* 

m o  = 


16  1  sin(^ix/2)sin(6ii/2) 

7T  K  +  1  <5? 


and  6x  is  given  by  (22b).  Subtracting  from  the  integrand  in  (104),  evaluating  some 

of  the  integrals  in  closed  form,  using  an  integration  cutoff  point  for  the  infinite 
integral  and  taking  the  limit  as  y  — >  0,  after  some  manipulations  can  be  written  as, 


(x,  0,  t )  =  («  +  1)  {  r«  (x,  t)  +  4fc)  (*,  t) 


miexp(7(t  -  x)/2)Ei(  -  A$(t  +  x))  j. 


(108) 


where  Ei()  is  the  exponential  integral  [10]  and, 


r22  (x>  t) 


+  - - i-o  +  — —  ?exp(7 (t  -  x)/2), 


(t  +  x)6  (t  +  xY  t  +  x 


(109a) 


fA% 

\x,t)=  Kft*(a,t,x)da 

JO 

poo  r  - 

+  /  t,x)  —  °°(a:,  f,x)  da 

L  J 

poo 

+  jA(hr^(a,t,x)da, 


(109b) 


t,  x)  =  K^\a,  t,  x)  -  [m*2a2  +  m\a  +  m0)exp(7(f  -  x)/2  -  (t  +  x)a), 


(110a) 


t,x)  =  {^  +  — |  H - f  lexp(7(f  -  x)/2  -  (t  +  x)a).  (110b) 

[  a2  a6  a1  J 


Second  integral  in  (109b)  will  be  neglected  in  numerical  computation  for  a  sufficiently 
large  value  of  A$.  Third  integral  in  (109b)  is  evaluated  in  closed  form.  The  expression 
used  in  the  evaluation  of  this  integral  is  given  in  Appendix  A.  k22  (x,  0,  t)  and  k$  (x,  0,  t) 
are  given  by  equations  (102)  and  (108)  respectively.  Adding  these  two  equations 


fc22  (x,  0,  t)  is  expressed  in  the  following  form, 
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k22(x,0,t)  =  exp(7x){  —  — — - f  h22s(x,t)  +  h22f(x,t) 

7 Xt  —  X 


(ill) 


where, 

h22s(x,t)  =  ^Y^r${x,t),  (H2a) 

h22f(x,t)  =  {k  +  1)  —  |  —  fell  In  ^221 1  £  —  ^1^  +  ^221  (X>  t)  t^222(X>^)|' 

+  ~~Y~{J^>^x't^ -miexp(7(i -x)/2)Ei(- A^  +  z))  j-  (112b) 

The  terms  in  (112)  are  given  explicitly  by  (103)  and  (109).  Note  that  the  first  term  in 
(1 1 1)  is  the  Cauchy  singularity  associated  with  a  crack  in  infinite  plane,  as  for  the  second 
term,  we  can  write, 


lim  h22s(x,t) 


if  1  2  £ _ 4£2 

7T  l  t  +  x  +  (t  +  x)2  (£  +  a:)3 


(113) 


This  term  becomes  singular  as  x  and  £  simultaneously  approach  zero,  and  is  the  standard 
expression  found  for  edge  cracks  in  homogeneous  materials,  (see  for  example,  equation 
(23a)  in  Dag  and  Erdogan  [8]).  Also,  note  that  if  the  medium  is  homogeneous  (i.e., 
7  =  0),  h22f  =  0  in  (111)  and  hnf  =  0  in  (83)  and  fcn(»,0,£)  =  k22(x,0,t),  but  if 
7^0,  this  equality  is  not  valid  and  consequently  for  graded  materials 
ku(x,0,t)  #  k22(x,0,t). 


3.4  k23(x,  y,  t ) 

Referring  to  (16c),  k22(x,  y,  £)  can  be  written  as, 

V >  *)  =  J  Mp>  x)exp(ip(y  -  t))dp, 


(114a) 
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(114b) 


023  (p,  x)  =  ^(*p  +  iV;sj)V’;(p)exp(sJx), 

j= 1 

where  Sj,  (j  =  1,2, 3, 4),  JV,-  and  ipj(p)  are  given  by  (13),  (15)  and  (19),  respectively. 
Changing  the  limits  of  integration  in  (114a)  and  taking  the  limit  as  y  — *  0,  and 
rearranging,  k23(x,  0,  t )  can  be  written  as 


k23(x,0,t)  =  6XP^— ^  {  K231(p,x)cos(pt)  -  K232(p,  x)sin 


(115) 


where. 


K23i(p,x)  =  023 (P,  x)  +  023 (  “  P,  x), 


K232(p,  X)  =  i(023(p,  -  023 (  -  P,  *))• 


(116a) 


(116b) 


In  order  to  extract  the  singular  terms  we  expand  K23i(p,  x)  and  K232(p,  x)  into  series  as 


p  — »  oo.  Following  asymptotic  expansions  are  obtained  by  using  MAPLE, 

■^231  (P> x)  =  ie\iP  +  eio  +  —  +  -y  H - f  — y  |exP(  ~  (7/2  +  p)x), 

l  P  P  P  ) 

K^2 (p,  z)  =  |e2iP  +  e20  +  ”  +  -Jf  4 - F  |exP(  —  (t/2  +  p)x), 

where  the  coefficients  are  functions  of  x.  The  leading  terms  are  given  by, 


en  —  — 


4:rjsin(83x/2) 


4sin(^3x/2) 


(117a) 


(117b) 


■5  C21 


(118a,b) 


and  83  by  (12b).  Subtracting  the  asymptotic  expansions  from  the  integrands  in  (115), 
using  integration  cutoff  points  for  the  infinite  integrals  and  evaluating  some  of  the 


integrals  in  closed  form,  after  some  manipulations  (115)  is  reduced  to, 
k23(x,0,t)  =  exp(7x)|/i23s(a:,t)  +  h23f(x,t)}, 


(119) 


h2zs{x,t)  =  ^{r23i(x,f)  +  r232(x,t)}, 

^23  l/(x,  f)  =  ^  |  ^231 0  +  ^232  + 


(120a) 


en exp (  -  jx/2) 


[r(0,  ( x  —  it)A23i)  +  r(0,  ( x  +  it)A2zi)] 


+  e2iexp(  —  7x/2)arctan(  —  £/x)  1, 


(120b) 


where  A2^i  is  an  integration  cutoff  point,  T(, )  is  the  incomplete  gamma  function,  [10], 


r231<1’!)  =  { (^fe'n  +  ?f?ei0}exp(  ” lx/2)’ 
r232<:M)  =  {j0keh  "  ^Tpe20}exp(_7x/2)’ 

M  231 

Jm(x,t)=  K2U(p,  x)cos(pt)dp 
Jo 

poo 

+  /  (^231  (p,x)  -  K^(p,x))cos(pt)dp 

J  A  231 
poo 

+  /  A1Up>*)cos(p*)dp, 

J  ^-231 

pA  232 

J2Z2{p,x)=  -  /  -^232 (P>  x)sin(pt)dp 

Jo 

poo 

-  /  (i^232(p,  x)  -  ^(p.  x))sin(pf)dp 

j  A2Z2 

poo 

-  /  r^2(p,x)sm(pt)dp, 

J  Am 

and  A 232  is  another  integration  cutoff  point.  The  remaining  terms  are  given  by 


(121a) 


(121b) 


(121c) 


(121d) 


K^(p,x)  =  KMp,x)  -  (e*np  +  eio)exp(  -  (7/2  +  p)x), 


(122a) 
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r&(p,  x)  =  { ^  ^  +  •• ••  +  ~§ }exp(  -  (7/2  +  p)x), 

K%32 (p,  x)  =  K232(p ,  x)  -  (e*2lp  +  e20  +  e2i/p)exp(  -  (7/2  +  p)x), 

r&M  =  {j  +  j  +  • -  - +  ^}-p(  -  (7/a  +  (>)*)• 


(122b) 

(122c) 

(123d) 


Second  integrals  in  (121c)  and  (12 Id)  will  be  neglected  in  numerical  computation  for 
sufficiently  large  values  of  A2 31  and  A2 32-  Third  integrals  in  (121c)  and  ( 1 2 1  d)  are 
evaluated  in  closed  form.  The  expressions  used  in  the  evaluation  of  these  integrals  are 
given  in  Appendix  A.  Also  note  that 


lim  h23s(x,t)  =  - 
>0  7T 


2  xt2 

"V+t2)2 


2  tx2  \ 

(x2  +  f2)2  J 


(124) 


This  term  becomes  singular  as  x  and  t  simultaneously  approach  zero,  and  is  the  standard 
expression  obtained  for  the  homogeneous  materials  (see,  for  example,  equation  (23c)  in 
Dag  and  Erdogan  [8]). 


3.5  k3i(x,y,  t ) 

We  first  express  k3i(x,  y,  t )  as  follows: 

hi{x,  y,  t )  =  k$(x,  y,  t )  +  k^\x,  y,  t ),  (125) 

where  k3^  and  k^  are  obtained  from  the  infinite  plane  and  half-plane  (x  >  0)  solutions 
respectively.  Referring  to  (25f),  ki^  (x,  y,  t)  can  be  expressed  as, 

1  r°° 

k$(x,  y,t)  =  -  y)exp(iu(x  -  t))du,  (126a) 

Tf  J  -00 

031  (w,  V )  =  ^PJ(w)nJexp(nJy) ,  (126b) 

3= 3 
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rij,  (. j  -  1, 2, 3, 4),  Pj  are  given  by  (23)  and  (29),  respectively.  Changing  the  limits  of 
integration,  (126a)  can  be  written  as 

kz\(x,y,t)  =  |  y)cos(ui(x  -  t))  +  K^2(u,y)sm(u(x  -  i))}dw  (127) 

where 

kSm  y)  =  *$(«,  y )  +  4i(  -  y),  (128a) 

-^312  (w>  y ) =  *  (^3?  (w>  y) _  03i  ( -  y))  •  (128b) 

In  order  to  extract  the  singular  terms  we  expand  and  K^2  into  series  as  u  —*  oo. 
Following  asymptotic  expansions  are  obtained  by  using  MAPLE 

•^311°° (w>  y)  =  { 9*nu  +  9io  +  —  +  H - F  “jj  |exp(  —  wy))  (129a) 

^  u)  u;  cj  j 

tfiuV, iO  =  {<>>  +  S2«  +  ^  +  ■  ■  ■  +  }exp(  -  w»),  (129b) 

where  the  leading  terms  are  given  in  Appendix  B.  Note  that  the  coefficients  of  the 
expansion  are  functions  of  y.  Subtracting  the  asymptotic  expansions  from  the  integrands 
in  (127),  using  integration  cutoff  points  for  the  infinite  integrals,  evaluating  some  of  the 
integrals  in  closed  form  and  taking  the  limit  as  x  —*■  0,  after  some  manipulations 
fcgi  (x,  y,  t)  is  reduced  to 

kl'!(0,y,t)  =  ~  jrji'i'iA  t)  +  JziiU),  t)  +  J-iMy,  t) 

+  [r(o,  (y  -  it)A^)  +  r(o,  (y  + 

+  <?2iarctan(  —  t/y) 1,  (130) 
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where  is  an  integration  cutoff  point,  T(, )  is  the  incomplete  gamma  function  [10], 


a),  y2-t2  *  ,  y  2ty  ^  1 

r3.  (»•  *>  =  +  JT¥gi° "  (7T??521 "  V2  +  t2 '■ 


L  (V,t)  =  f 
JO 


Kziiiw,  y)cos(ut)du 


+  f  ;)  {^(w,  y)  -  (u,  y)  }cos(ut)do 


(131a) 


.  r^(u,y)cos(ut)du, 


(131b) 


Jm(v>t)=  ~  /  Kg£(w,y)sm(ut)dw 
Jo 


f  w  {^312(^ 2/)  -  ^312° 0*>,  y)}sin(ut)du 


r3^2  (a;,  y)sin(utf  )du;, 


(131c) 


Ai*, I  is  another  integration  cutoff  point  and  the  remaining  terms  are  given  by 


#311  (w»  V )  =  ^3ii(w»  2/)  -  (0nw  +  0io)exp(  -  wy), 

4;f>> »)  =  { §r + ^ + •  ■ ■  • + }exp(  -  ^). 


(132a) 


(132b) 


(w,  y)  =  if 312(^3  2/)  -  (02iw  +  520  +  W^)exp(  -  uy), 


r& $>.»)  =  { 


922  923 

u2  uz 


+  ^77  lexP(  -  w2/), 


(132c) 


(132d) 


Second  integrals  in  (131b)  and  (131c)  will  be  neglected  in  numerical  computation  for 
sufficiently  large  values  of  and  A^2  ■  Third  integrals  in  (131b)  and  (131c)  are 
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evaluated  in  closed  form.  The  expressions  used  in  the  evaluation  of  these  integrals  are 
given  in  Appendix  A.  Referring  to  (35f)  and  (41a),  k^(x,  y,  t )  can  be  written  as  follows: 


poo 

k£\x,  y,  t)  =  2  /  i,  x)sin (ay)da, 

Jo 


(133) 


where 


K^{a,t,x)  = 


-  a^B*(a,t)exp(pjx  +  (7/2  -  Ai)f), 


(134) 


pj,  (j  =  1,2, 3, 4)  and  B*(a,t)  are  given  by  (34)  and  (41),  respectively.  In  order  to 
extract  the  singular  terms  we  expand  K^\a,t,0)  into  a  series  as  a  tends  to  infinity  as 


follows: 


K^)oo(a,t,  0)  =  |iia  +  i0  +  ^  +  ^  +  ---  +  ^f|exp((7/2-a)t)  (135) 

where  the  coefficients  of  the  expansion  are  functions  of  t.  The  leading  term  is  given  by 


2  «;sin(<5if/2) 

7T  (K  +  l)6i 


(136) 


61  by  (22b).  Subtracting  K^}o°  from  the  integrand  in  (133),  evaluating  some  of  the 
integrals  in  closed  form,  using  an  integration  cutoff  point  for  the  infinite  integral  and 
taking  the  limit  as  x  — *  0,  after  some  manipulations  can  be  written  as 


k?M  y,  t)  =  2|r^}(y,  t)  +  iiexp(7f/2)arctan(y/f)  +  J$\y,  i)}, 


(137) 


(138a) 
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sin(m/)da: 


(?»  t)=  K31*  t)sin(ay)da 

Jo 

rOC 

+  Mtf'fW.o)-  ^“(a.t.o) 

^31 
poo 

+  l^r^°°(a,t)sm(ay)da.  (138b) 

JA3i 

In  (138)  Afi  is  an  integration  cutoff  point  and  the  remaining  terms  are  given  by 
K^*(a,t)  =K£\a,t,0)-  (i[a  +  i0  +  ii/a)exp((7/2  —  a)*),  (139a) 

^(M)  =  {I  +  |  +  -  +  |}«P((7/2  -  a)t).  (139b) 

Note  that  second  integral  in  (138b)  will  be  neglected  in  numerical  computation  for  a 
sufficiently  large  value  of  and  third  integral  is  evaluated  in  closed  form.  The 
expression  used  in  the  evaluation  of  this  integral  is  given  in  Appendix  A.  *4?  (0,2/,*)  and 
k^(0,y,t)  are  given  by  equations  (130)  and  (137),  respectively.  Adding  these  two 
equations,  kzi  (0,  y,  t)  is  expressed  in  the  following  form, 

hi(0, y,t)  =  hu(y,t)  +  h31f(y,t),  (14°) 


where, 

hu{y,t)  =  -rzi(y,t)  +  2r{^(y,t), 

7T 


(141a) 


9n 

2 


+ 


r(o,  {y  -  it)A\ jjj)  +  r(o,  (y  +  i*)4n)  j  +  02iarctan(  -  t/y) 
2|4i)(y>0  +  *iexp(7*/2)arctan(y/*)|. 


(141b) 
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The  terms  in  (141)  are  given  explicitly  by  (131)  and  (138).  Also,  note  that, 


lim  h3u(y,t) 
M-* 0 


4  _fy_ 
n  (y2  +  t2)2' 


(142) 


This  term  becomes  singular  as  y  and  t  simultaneously  tend  to  zero,  and  is  the  standard 
expression  found  for  homogeneous  materials,  (see,  for  example,  equation  (23d)  in  Dag 
and  Erdogan  [8]). 


3.6  k32(x,y,  t) 


We  first  express  k32(x,  y,  t )  as  follows, 

k32  (x,y,t)  =  k^(x,y,t)  +  k$(x,y,t),  (143) 

where  k^  and  k^  are  obtained  from  the  infinite  plane  and  half-plane  (x  >  0)  solutions 
respectively.  Referring  to  (45 f)  k22  (x,  y ,  t)  can  be  expressed  as 
i  r°° 

k$(x,  y,t)  =  -  y)exp(iu(x  -  t))du,  (144a) 

J -ex) 

032  (w>  y )  =  y^zi{u)nitxv{njy),  (144b) 

i= 3 

where  rij,  ( j  =  1, 2, 3, 4)  is  given  by  (23)  and  Zj  is  given  by  (48).  Changing  the  limits  of 
integration  (144a)  can  be  written  as 

k$(x,  y ,  i)  =  ^  jf  {  Kz2i(^ y)cos(uj(x  -  t))  +  K$2(u,  t/)sin(u>(a;  -  t))}du,  (145) 

where 

^321  (w>  y)  =  0S  (w>  y)  +  0S  ( -  <*>>  y)> 


(146a) 


#312 (w>  V )  =  i (^32  (w>  ?/)  “  032  (  ~  2/))  •  (146b) 

In  order  to  extract  the  singular  terms  we  expand  K^i  and  into  series  as  u  — *  oo. 
Following  asymptotic  expansions  are  obtained  by  using  MAPLE: 

T(U’V)  =  l^iiW  +  iio  4-  ~  +  “I  H - h  |exp(  -  uy),  (147a) 

y)  =  +  *20  +  ^  ^  ■ ■  +  jf}exp(  -  uy),  (147b) 


where  the  leading  terms  are  given  in  Appendix  B.  Note  that  the  coefficients  of  the 
expansion  are  functions  of  y.  Subtracting  the  asymptotic  expansions  from  the  integrands 
in  (145),  using  integration  cutoff  points  for  the  infinite  integrals,  evaluating  some 
integrals  in  closed  form  and  taking  the  limit  as  x  — *  0,  after  some  manipulations 
^32  (x,  y,  t )  is  reduced  to: 


*32  (0,  y,  t ) 


+ 


‘ii 


r(o,  (»  -  +  r(o,  (y  + 


+  Z2iarctan(  -  t/y) 


5 


(148) 


where  A^\x 

r&fat)  = 


is  an  integration  cutoff  point,  T(, )  is  the  incomplete  gamma  function, 
y2-t2  „  ,  y  7  2 ty  ,*  t  , 

(y2 +^4" +  3,2  +  12  (3,2 +(2)2 '2'  j,2  +  f2  '=0- 


(149a) 
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A322  is  another  integration  cutoff  point  and  the  remaining  terms  are  given  by 


#f> (w,  y)  =  #321  (^>  y)  -  (l*nu  +  iio)exp(  -  coy), 

r3?r(".  v)  =  {  £5  +  +  •  •  •  +  jf  }exP(  - 


(149b) 


(149c) 

(150a) 

(150b) 


#322  (w,  y)  =  #322  (W>  1/)  -  (z21w  +  *20  +  W^)eXP(  ~  Wy), 


#322  (w>  y) 


u;2  ur  a;12  I 


(150c) 

(150d) 


Second  integrals  in  (149b)  and  (149c)  will  be  neglected  in  numerical  computation  for 
sufficiently  large  values  of  A§ \  and  A§2.  Third  integrals  in  (149b)  and  (149c)  are 
evaluated  in  closed  form.  The  expressions  used  in  the  evaluation  of  these  integrals  are 
given  in  Appendix  A.  Referring  to  (54f)  and  (61a),  k$(x,  y,  t )  can  be  written  as  follows: 

POO 

k${x,y,t)=2  #32 ^  x)cos(ay)da,  (151) 

Jo 
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where 


K${a,t,x)  =  aJ2G*(a,t)txTp(pjX  +  (7/2  -  \{)t), 
j= 3 


(152) 


tj,  (j  =  1, 2, 3, 4)  is  given  by  (53)  and  G*(a,  t)  is  given  by  (60).  In  order  to  extract  the 
singular  terms  we  expand  K^{cn,  t,  0)  into  a  series  as  a  tends  to  infinity  as  follows 

=  Ua  +no  +  ^  +  2f  +  •••  +  ^}exp((7/2  -  «)<),  (153) 

where  the  coefficients  of  the  expansion  are  functions  of  t.  The  leading  term  is  given  by 


ni  =  ~  - 


2  /csin(#it/2) 


(154) 

7T  (k  +  1)^1 

and  Si  by  (22b).  Subtracting  k£)c°  from  the  integrand  in  (151),  evaluating  some  integrals 
in  closed  form,  using  an  integration  cutoff  point  for  the  infinite  integral  and  taking  the 
limit  as  x  — »  0,  after  some  manipulations  can  be  written  as 

k£\0,y,t)  =  2  {r£\y,t)  +  jg\y,t)  + 

r(o,  (t  +  iy)A\ +  r(o,  (t  -  exp(7*/2)  j,  (155) 


ni 

2 


where, 

r« ' M  =  {  +  ^n}exP^/2)'  (156a) 

Aw 

=  /  32  K{^)*(a,t)cos{ay)da 
Jo 

+  r  [ffg  W,  0)  -  (o,  t,  0)]  cos(ay)da 

poo 

+  /  r3(2)oo(a,f)cos(ay)dQ.  (156b) 
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A%2  is  an  integration  cutoff  point,  and  the  remaining  terms  in  (156b)  are  given  by 
K^*(a,t)  =  K$\a,t,  0)  -  (nja  +  n0)ex p((7/2  -  a)t),  (157a) 


4‘ >“(<*,4)  =  {-!+-,+-  +  ^  W(t/2  -  «)t).  (157b) 

[  or  cr  cr  J 

Note  that  second  integral  in  (156b)  will  be  neglected  in  numerical  computation  for  a 
sufficiently  large  value  of  and  third  integral  is  evaluated  in  closed  form.  The 
expression  used  in  the  evaluation  of  this  integral  is  given  in  Appendix  A.  (0,  y,  t)  and 
k$(0 ,y,t)  are  given  by  equations  (148)  and  (155),  respectively.  Adding  these  two 
equations,  fc32(0,  y,  t)  maybe  expressed  in  the  following  form, 

MO  ,y,t)  =  h32s(y,t)  +  hZ2f(y,t),  (158) 

where 

h32s(y,t)  =  -rg(y,  t)  +  2r{$(y,t),  (159a) 

7 r 

h32  =  -  |«42i  (y>  t)  +  J^iVi  t)  + 

+  y  T(o,(y-  it)A$ i)  +  r(o,  (y  +  j  +  Z2iarctan(  -  t/y) J 

+  2| 4fc)(»,t)  +  y  r(o,(t  +  iy)Ag})  +r(0,(t-zy)Ag))]exp(7t/2)}.  (159b) 


The  terms  in  (159)  are  given  explicitly  by  (149)  and  (156).  Also,  note  that, 


lim  h32s{y,t) 

(y,*)— "0 


4  i3 

7T  (y2  _|_  *2)2  ' 


(160) 
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This  term  becomes  singular  as  y  and  t  simultaneously  tend  to  zero,  and  it  is  the  standard 
expression  found  for  homogeneous  materials,  (see,  for  example,  equation  (23e)  in  Dag 
and  Erdogan  [8]). 

3.7  k33(x,y,  t ) 


Referring  to  (16d),  fc33(x,  y,  t)  can  be  written  as 
2  1  f°° 

k33(x,  V >  *)  =  —T  ~  x)exp(ip(y  -  t))dp, 

K  +  1  7T  J-oo 


(161a) 


Mp,x)  = 

j= 1 


(161b) 


where,  (j  =  1,2, 3, 4)  and  0j(p)  are  given  by  (13)  and  (19),  respectively.  Changing 
the  limits  of  integration  in  (161a),  £33  can  be  written  as 

k33 (x,y,t)  =  —  [#33i(p,x)cos(p(y  -  t ))  +  #332(p,x)sin(p(y  -  t ))  dp, 


(162) 


where. 


#331  (p,  ®)  =  033  (p,  +  033  (  -  P,  *), 


(163a) 


#332  (Pj  ®)  =  *(033(P,  *)  -  033 (  ~  P,  *))• 


(163b) 


In  order  to  extract  the  singular  terms  we  expand  #331  and  #332  into  series  as  p  — *  00. 


Following  asymptotic  expansions  are  obtained  using  MAPLE: 


#331  (p,x)  = 


K&fax) 


.  Cn  C12 

C10  H - 1 — ^  + 

P  P 2 


.  c2i  c22 

C20  H - F  -r  + 

p  p2 


+  ^f}exP(“Px)’ 


(164a) 


(164b) 


49 


where  the  leading  terms  are  given  by 


Cio  ="??(«- 1)>  c20  =  ^(k  +  1).  (165a, b) 

Subtracting  the  asymptotic  expansions  from  the  integrands  in  (162),  using  integration 
cutoff  points  for  the  infinite  integrals,  evaluating  some  integrals  in  closed  form  and  taking 
the  limit  as  x  — »  0,  after  some  manipulations  (162)  is  reduced  to: 

kx(0,y,t)  =  (  -  !-!-  +  (166a) 

[  nt  —  y  k  + 1  J 

2  f  7T 

h33f(y,t)  =  -T-— r  <  -  culn(^33i|i  -  2/1)  ~  c2iTsign(t  -  y) 

7 T{K  +  L)  l  l 

+  7331(2/5 1)  +  J332(y>  t) | ,  (166b) 


where  _A331  is  an  integration  cutoff  point,  <5()  is  the  Dirac  delta  function. 


/  ^331  r 

=  J  [#331 0, 0)  -  Ciojcos(p(y  -  t))dp 
+  JA  [#331  (p,  0)  -  #3“i(p,0)]cos(p(y  -  t))dp 
+  f  [#331 (P>°)  -  cio  -  cn/p  cos (p(y  -  t))dp 

J  ^331 

(  3/1  cos(o:)  —  1  1 

—  cn<  70  +  /  - - - da}, 


^332(2/)  t)  =  J  [#332 (p>  0)  -  c2o  -  C2i/pj  sin (p(y  —  t))dp 

poo  r 

+  /  #332 (P)  0)  -  #332 (P)  °)j  sin(p(2/  -  0)dP 

^  AzZ2  1 

poo  r  -j 

+  /  [#^2(p)°) -c20-C2i/pjsin(p(y-t))dp, 

v  A332  L 


(167a) 


(167b) 
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A332  is  another  integration  cutoff  point  and  70  is  the  Euler  number,  [10].  Second  integrals 
in  (167a)  and  (167b)  will  be  neglected  in  numerical  computation  for  sufficiently  large 
values  of  A331  and  Aw  Third  integrals  in  (167)  are  evaluated  in  closed  form.  The 
expressions  used  in  the  evaluation  of  these  integrals  are  given  in  Appendix  A.  Also,  note 
that  first  two  terms  in  (166a)  are  the  singular  terms  and  they  are  the  standard  expressions 
obtained  for  homogeneous  materials,  (see,  for  example,  equation  (22c)  in  Dag  and 
Erdogan  [8]). 

The  contact  problem  for  decreasing  stiffness  (7  <  0) 


16t]{k2(1  +  sign(7))  -  /c(7  +  9sign(7))  +  10  +  16sign(7)} 
74(k  +  1)2(1  +  sign(7)) 


(170b) 


__  2(/c  -  1) 
01  70  +  1)  ’ 


(170c,d) 


8(fc  -  3)(k(1  +  sign(7))  -  4) 

73(l  +  /c)2(l  +  sign(7)) 

Observing  that  sign(7)  =  1  for  7  >  0,  sign(7)  =  —  1  for  7  <  0  from  (169)  and  (170)  it 
is  seen  that  Kzzi(p,  0)  and  K332(p,  0)  are  well  behaved  near  p  =  0,  for  7  >  0.  However, 
for  7  <  0  coefficients  64  and  a3  (and  possibly  that  of  higher  powers  of  p)  become 
unbounded  and  as  a  result  £33  expressed  by  (166a)  also  becomes  unbounded. 
Consequently,  it  is  seen  that  for  a  graded  half  plane  with  an  exponentially  decaying 
stiffness  the  contact  problem  is  not  a  well-posed  problem.  Physically,  the  problem  that  is 
analogous  to  7  <  0  case  is  a  homogeneous  strip  of  finite  thickness  under  an  unbalanced 
transverse  load  P  (in  thickness  direction)  which  has  no  solution  (see  Ratwani  and 
Erdogan  [9],  for  explanation).  Thus  for  graded  half-planes  with  or  without  a  crack  if 
7  <  0  the  contact  problem  has  no  solution. 

4.  Singular  behavior  of  the  solution 

The  integral  equations  of  the  problem  are  given  by  (65)  and  the  kernels  of  the 
equations  are  derived  in  Section  3.  The  asymptotic  behaviors  of  the  integrands  are  also 
examined  and  singular  terms  are  extracted.  Using  the  expressions  given  in  section  3, 
integral  equations  given  by  (65)  can  be  written  as  follows: 

(Tyy(x,  0)exp(  —  7X)  =  J  + hUs{x,t)  +  hnf(x,t)^fi(t)dt  + 

+  [  hi3s(x,t)  +  hizf(x,t)  fz(t)dt  +  E0e0  =  0,  0  <  x  <  d,  (171a) 

J  a  L  J 

axy(x,  0)exp(  -  7«)  =  J  +  ^(M)  +  h22/(x,t)]/2(t)dt  + 

+  f  h23S{x,t)  +  h2zf(x,t)  fz(t)dt  =  0, 

J  a  J 


0  <  x  <  d, 


(171b) 


Y~ J^(0, y)  =  [^31  s{y,t)  +  hif{y,t)^fi(t)dt  + 

rd  r  ] 

+  /  hz2s  (y,  t )  +  hz2}{y,  t )  h{t)dt  — 

Jo  L  J 


-£r>>+. 


nt-y 


+  h33 f(y,t)  f3(t)dt,  a<y<b , 


(171c) 


where  expressions  for  hijS(x*,t)  and  hijf(x*,t )  (x*  =  x  for  i  =  1,2  and  x*  —  y  for 
i  =  3)  are  given  in  Section  3.  hijs(x*,t)  are  the  generalized  Cauchy  kernels  (of  the  order 
1/f)  that  become  unbounded  as  the  arguments  x*  and  t  tend  to  the  end  point 
simultaneously.  hijf(x*,t)  are  bounded  Fredholm  kernels.  The  singular  terms  are  found 
to  be: 

,  ,  1(1  2 1  U2  \ 

lim  hus(x,t)  =  lim  h22S{x,t)  =  —<  — - 1-  - - —  — - 3 

0  v  '  (x,t)—*o  Tr{t  +  x  (: t  +  x y  (t  +  xy  J 


(x,t)—>0 


0  <  (t,  x)  <  d, 


1  f  2  xt2  2  i3  1 

}■  “  <  ‘  <  M  <  *  «*, 


(172a) 


(172b) 


,  If  2  xt2  2tx2  1 

lim  h23S(x,t)  =  -<  77— - —5  -  7— - rr? 

(*,*)— 0  v  ir\'(t2  +  x2)2  (t2  +  x2)2  J 

hm  h31s{y,t)  =  ~  ~  .2, 

(v.*)-0  7T  (<2  +  y2) 

1  4t3 

lim  h32s{y,t )  =  -  2r2~, 

(y.t)-*0  7T  (j?  +  yi) 


a  <t  <b,  0  <  x  <  d,  (172c) 


0  <t<d,  a  <  y  <b,  (172d) 


0  <  t  <  d,  a  <  y  <b,  (172e) 


Note  that  the  singular  terms  in  the  integral  equations,  i.e.,  the  Cauchy  singularities  and 
generalized  Cauchy  kernels  given  by  (172)  are  also  obtained  for  the  crack/contact 
problem  in  a  homogeneous  half-plane.  If  we  compare  (172)  and  equations  (23a, e)  which 
are  given  in  [8],  we  observe  that  the  singular  terms  are  identical  except  for  the  sign 
changes  for  some  terms,  which  are  due  to  the  different  coordinate  axes  used  in  [8]  and  in 


this  report.  It  may  then  be  concluded  that  the  singular  behavior  of  the  unknown  functions 
for  the  graded  and  homogeneous  materials  are  identical.  Hence,  for  the  graded  materials, 
the  singular  behavior  of  the  unknown  functions  is  independent  of  the  material 
nonhomogeneity  constants  'y  and  /iQ  and  depend  on  the  friction  coefficient  ?/  and  the 
surface  value  of  the  Poisson's  ratio  (through  elastic  constant  /c)  only.  The  details  of  the 
function-theoretic  analysis  to  determine  the  singular  behavior  of  the  unknowns  are  given 
in  [8].  Here  we  summarize  the  results,  for  the  two  cases  a  >  0  and  a  —  0. 

a  >  0 

In  this  case  the  kernels  (172b-e)  are  bounded  in  their  corresponding  closed  intervals  and 
would  not  contribute  to  the  singularities  of  the  functions  fa,  fa  and  fa.  We  express  fi  as, 


fa(x)=xe'(d-x)x'F1(x), 

0  <  x  <  d, 

(173a) 

fa(x)  =  xd>(d-x)X2F2(x), 

0  <  x  <  d, 

(173b) 

Mv)  =  (y-  aT(b  -  y)0F3(y), 

a  <y  <b. 

(173c) 

The  function-theoretic  analysis  to 

determine  the  exponents  is 

described  in  [8],  and 

following  equations  are  obtained, 

01  =0,  02  =  0, 

(174a) 

cot(7rAi)  =  0,  cot(7rA2)  =  0, 

LO" 

O 

1 

II 

CM 

II 

(174b) 

^ 

cot(7ru;)  =  t] - -,  cot(7 r/3)  = 

AC  ~h  1 

AC  —  1 

=  \  +  i’ 

(174c, d) 

where  acceptable  roots  are  Ai  =  —  0.5,  A2  =  —  0.5,  3ft(u>)  <  0  if  a  is  known  and  is  a 
sharp  comer,  3ft  (gj)  >  0,  if  a  is  unknown  and  is  a  point  of  smooth  contact.  Similarly 
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&(/?)  <  0  if  b  is  a  known  sharp  comer  and  &(/?)  >  0  if  b  is  unknown  and  the  contact  is 
smooth. 

a  =  0 

In  this  case  all  kernels  hijs(x*,t )  become  unbounded  as  x*  — ►  0  and  t  -»  0 
simultaneously  and  contribute  to  the  singularity  of  the  unknown  functions.  Again,  we 


express  the  unknown  functions  as  follows, 

h{x)  =  xa(d  -  x)AlGiOr),  0  <x<d,  (175a) 

f2(x)  =  xa(d-x)x'G2(x),  0  <  x  <  d,  (175b) 

fz(y)  =  ya(b-y)0Gs(x),  0  <y<b,  (175c) 

The  function-theoretic  analysis  carried  out  in  [8]  shows  that 

cot(7rAi)  =  0,  cot(7rA2)  =  0,  (Ai  =  A2  =  —  0.5),  (176a, b) 

cot(7 r/3)  =  —  77— — (176c) 


As  shown  in  [8]  equation  (56)  and  (57),  the  eigenvalue  a  and  the  expressions  relating 
G2(0)  and  G 3(0)  are  given  by 

2a2  +  4a  +  1  -  cosset) 

( k  +  l)sin2(7ra) 

x  (?7(4o;2  +  10a  +  5  +  (k  -  l)cos(7ra)  +  k( 2a  +  3))  +  (k  +  l)sin(7ra))  =  0,  (177a) 


Gi(0)y/d 


rjja  +  2)cos(7ro:/2)  +  (a  +  l)sin(?ra/2)  }  Q 
2 a2  +  4a  +  1  -  cos(7ra)  /  3 


(177b) 


(177c) 
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From  the  symmetrically  loaded  half  plane  (y  >  0)  analogy  it  is  known  that  a  =  0  for 
r/  =  0  and  at  x  =  0,  y  =  0  the  stress  state  is  bounded.  In  this  case  function-theoretic 
analysis  shows  that  in  order  not  to  have  a  logarithmic  singularity  in  the  integral  equations 
following  condition  must  be  satisfied, 

G2(0)Vd=^G3(0)b13.  (178) 

5.  Numerical  solution  of  the  integral  equations 

In  this  section,  we  will  develop  a  numerical  solution  method  for  the  case  of  a  flat 
stamp  as  shown  in  Figure  2.  We  first  normalize  the  intervals  in  (171)  by  defining 


P 


Figure  2:  The  geometry  of  the  crack/contact  problem  for  a  flat  stamp. 


d  d 
t=2r+2’ 


i  = 


b  —  a 
2 


r  + 


b  4-  a 
2 


(179a, b,c) 


in  integrals  involving  /2(f)  and  /3(f),  respectively.  Then  we  define  the  normalized 
unknowns  of  the  problem  as  follows: 
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(180a) 


Mr) 


Mr) 


Mr)  = 


h  (r + f) 

PKb-a)  ’ 

h(tr+t) 

P/(b-a)  ’ 

,  (b  —  a  b  +  a 

h{—r+^r 

P/(b  -  a) 


-  1  <  r  <  1, 


—  1  <  r  <  1, 


1  <  r  <  1. 


(180b) 


(180c) 


The  intervals  (0,  d)  and  (a,  6)  are  also  normalized  by  defining, 


d  d 
x~2Sl+2’ 


x=2S2  +  2’ 


b  —  a  b  +  a 
2  53  2  ’ 


for  eqn.  (171a), 


for  eqn.  (171b), 


for  eqn.  (171c). 


(181a) 


(181b) 


(181c) 


Using  (180)  and  (181),  integral  equations  (171)  and  equilibrium  condition  (7)  can  be 


written  as 


1  f  Mr)  dr+  f  jjn{Sur)MT)dr  +  /  Hn(si,r)Mr)d7 

J-i  r  —  S\  J- 1 


- Eq€q -  q 

P/(b  -  a)  ’ 


-  1  <  Si  <  1, 


(182a) 


-  f  — —dr+  [  H22{s2,r)Mr)dr+  [  #23(s2,  r)Mr)dr  = 
7T  J-\T  —  S2  J- 1  7-1 


1  <  S2  <  1, 


(182b) 


"1  /*1  K  —  1 

Hn{si,r)Mr)dr  +  /  #32(53,  r)<j>2(r)dr  ~  r\— rr MM  + 

-1  7-1  k  +  1 

+  J  H3z(sz,r)Mr)dr  ~  0,  -  1  <  s3  <  1, 


(182c) 


J  H(s3)dsz  =  —  2, 


(182d) 


where  Hij(si,r )  are  given  in  Appendix  C.  Unknown  function  <fri(r),  (i  =  1,2,3)  can 
now  be  expressed  in  the  following  form, 


71=0 


71=0 


Hr)  = 


Wi (r)  =  (1  -  r)_1/2(  1  +  r)ai, 

(183a) 

w2(r)  =  (1  -  r)_1/2(l  +  r)a\ 

(183b) 

w2(r)  =  (1  -  r)P(l  +  r)“2, 

(183c) 

71=0 


where  for  a  =  0,  a\  —  a.2  —  a  and  for  a  >  0  ol\  =  0,  0.2  —  w.  Substituting  (183c)  into 
(182d),  A30  is  obtained  as 

2^+Q2+ir(/3+1)r(a2  +  1) 


^30  =  —  2/00) 


00  = 


r(/3  +  0:2  +  2) 


(184a,b) 


Now,  substituting  (183)  into  (182),  regularizing  the  singular  parts  of  the  equations  using 
the  expressions  given  in  Appendix  D  in  [8]  and  truncating  the  infinite  series  at  N, 
following  system  of  linear  algebraic  equations  is  obtained: 

N  N  ^  ^ 

Tmnn(s1)Aln  +  ^^m13n(si)A3n  =  -  — —  r  -  mi3o(si)A3o, 

n= 0  n=l  ^  '  '°  a > 

-  1  <  si  <  1,  (185a) 


71=0 


N  N 

^^^2271(52)^.277  +  J^jTn23n(S2)Azn  =  ~  W23o(s2)^-30!  —  1  <  S2  <  1, 

71=1 


(185b) 


N  N  N 

^^^3171(53)^171  +  "^^^3271(53)^271  +  ^^33?! (53)^371  ”  ”  ^330(53)^-305 

71=0  71=0  71=1 

-  1  <  s3  <  1.  (185c) 
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The  expressions  for  mijn{si),  (i,  j  =  1, 2, 3)  are  given  in  Appendix  C.  Note  that  if  a  =  0, 
compatibility  conditions  expressed  by  (177)  and  (178)  must  also  be  considered. 
Substituting  (183)  in  (177)  and  (178)  we  obtain, 


T)  >  0 

- 1)  -  h(4iY°2Hi,2tA^A)(  -  d = 

71=0  '  '  »=1 

=  f1(a)(D  2/3+1'/2A30, 

\a  ) 

jrA2nPt'M(  -  1)  -  A(a)(i)  -  1)  = 

n= 0  '  '  n=l 

=  h  2/J+1/,2A.3o, 

77(0:  +  2)cos(7ra/2)  +  (a  +  l)sin(7ra/2) 
f1'01'  ~  2a2  +  4a  +  1  -  cos(7rai) 

rj(a  +  l)sin(7ro:/2)  —  o;cos(7rQ:/2) 
j2{oc)  —  2a2  +  4a  +  1  —  cos(7tq;/2) 


(186a) 

(186b) 

(186c) 

(186d) 


77  =  0 

f>„P,<-1A0)(  -  1)  -  ')(  _  D  =  (187) 

n=0  n= 1 


Equations  (185)  can  be  solved  using  the  collocation  technique.  For  a  >  0,  the  number  of 
unknowns  is  (3 N  +  2).  Roots  of  the  Chebyshev  polynomials  are  used  as  the  collocation 
points  as  follows: 


Su  =  cos 


tt(2*-1)\ 
2(N  +  1)J' 


i  =  1, . . . ,  N  +  1, 


(188a) 


(  w(2t  —  1)  \ 

U^  +  l)/ 


S2j  =  cos 


i  =  1, . . . ,  N  +  1, 


(188b) 


—  COS 


tt(2 i  -  1) 


i  = 


(188c) 


In  the  numerical  solution  for  a  —  0,  the  equations  (186)  and  (187)  are  also  considered. 
After  solving  equations  (185)  for  Ain,  ( i  =  1,2,3)  the  contact  stresses  axx(0 ,y)  and 
<7Iy(0,  y)  and  stress  intensity  factors  at  the  crack  tip  (d,  0)  may  be  evaluated  by  using  the 
results.  The  stress  intensity  factors  are  defined  by  and  calculated  from 


kj  —  lim  y/2(x  —  d)ayy(x,  0)  = 

x^d+0 


=  -  lim  'l  y/: 2 (d  -  x)  ^  (v(x,  0+)  -  v(x,  0  )), 

x->d-0  K  +  1  O 


(189a) 


kn  —  lim  \/2(x  -  d)axy(x ,  0)  = 

x-*d+ 0 


=  -  lim  -  -  x)-^-(u(x,  0+)  —  u(x,  0  )). 

x— 0  K  "h  1  OX 


(189b) 


Using  (189),  the  normalized  stress  intensity  factors  and  the  normal  component  of  the 
contact  stress  may  be  expressed  as 


r  0  —  a*—' 


(190a) 


(190b) 


(  b  —  a  b  +  a 
a“(0'  _2_S3  +  “ 
P/(b  -  a) 


=  (i  -  + s3rj2A^3m>^- 


(191) 
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6.  Results  and  discussion 


The  calculated  results  in  this  report  consist  of  the  normal  and  in-plane  components  of 
the  stresses  on  the  surface  of  the  graded  half  plane  in  the  absence  of  a  crack,  contact 
stresses  (crxx(0,  y),  axy{0 ,y),  a  <  y  <  b)  and  the  stress  intensity  factors  (fc/,  kn).  First 
we  give  some  results,  showing  the  surface  stresses  in  a  graded  medium  in  the  absence  of  a 
crack  and  loaded  by  a  sliding  flat  stamp.  These  results  are  shown  in  Figures  3-8.  Contact 
stresses  in  the  absence  of  a  crack  are  calculated  by  solving  the  integral  equation  given  by 
(168)  and  also  considering  the  equilibrium  condition  given  by  (7).  As  shown  in  section 
3.7,  the  contact  problem  with  or  without  a  crack  has  no  solution  for  7  <  0.  Hence,  in 
Figures  3-8  results  are  given  for  positive  values  of  the  nonhomogeneity  parameter  7. 
Figures  3  and  4  show  that,  for  rj  =  0  both  axx(0,  y )  and  ayy( 0,  y)  are  symmetric  and  they 
have  square-root  singularities  at  y  =  a  and  y  —  b.  ayy( 0,  y)  vanishes  outside  the  contact 
area  for  7  =  0,  and  as  7  increases  it  becomes  tensile  at  both  ends  of  the  contact  region. 
Figures  5  and  6  show  the  results  for  77  =  0.4.  It  is  seen  that  there  is  a  greater  stress 
concentration  near  the  trailing  end  of  the  stamp  ,y  =  a  and  |<u|  >  |/3|,  u  and  /3  being  the 
singularities  at  y  —  0,  and  y  —  b  respectively.  The  important  conclusion  one  may  draw 
from  Figure  6  is  that  at  the  trailing  end  of  the  stamp  the  in-plane  component  of  the  stress 
<jyy(0,y)  is  unbounded,  tensile  and  discontinuous  and  has  a  singularity  of  the  order 
(a  -  y)w,  where  -  u>  1/2.  This  implies  that  y  =  a  is  a  likely  location  of  surface  crack 
initiation.  Similar  results  are  also  shown  in  Figures  7  and  8  for  77  =  0.8. 

In  Figures  9-16,  stress  intensity  factors  fcj  and  kn  are  shown  as  functions  of  the 
relative  stamp  size  b/d  for  a  =  0,  v  =  0.25  and  for  various  values  of  7  and  77.  In  these 
figures  the  stress  intensity  factors  are  normalized  with  respect  to  P/yfd  also  the 
nonhomogeneity  parameter  is  used  in  normalized  form  7 d.  The  circles  in  these  figures  are 
the  results  obtained  from  the  solution  of  the  homogeneous  half-plane  problem  as 
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described  in  Dag  and  Erdogan  [8].  It  is  shown  in  Figures  9  and  10  and  in  all  the  results 
presented  in  this  study  that,  the  results  obtained  from  the  solution  of  the  graded  half  plane 
problem  by  letting  7 d  =  0.0001  and  those  obtained  from  the  solution  of  the 
homogeneous  half  plane  problem  are  in  very  good  agreement.  That  is,  for  all  intents  and 
purposes  these  two  sets  of  results  are  identical.  Figures  9  and  10  show  that  for  a  =  0  and 
77  =  0,  i.e.,  for  the  case  of  normal  indentation,  mode  I  stress  intensity  factors  are  negative 
and  mode  II  stress  intensity  factors  are  positive  for  all  values  of  the  nonhomogeneity 
constant  7 d.  Mode  I  stress  intensity  factors  increase  and  mode  II  stress  intensity  factors 
decrease  as  7 d  increases.  Since  kj  is  less  than  zero  crack  closure  occurs  and  there  is 
contact  between  the  crack  faces.  But  the  results  can  still  be  applicable  and  useful  in 
superposition  with  an  uncoupled  solution  resulting,  for  example,  from  remote  strain 
loading  eyy(x,  Too),  [11],  [12],  provided  the  resultant  kj  is  positive.  Otherwise,  the 
problem  needs  to  be  formulated  by  taking  into  account  the  crack  closure  and  determining 
the  closure  distance  from  the  condition  of  kj  =  0.  Figures  11  and  12  show  the  modes  I 
and  II  stress  intensity  factors  for  77  =  0.2  and  for  a  =  0.  It  is  seen  that,  mode  I  stress 
intensity  factors  increase  as  the  friction  coefficient  increases  but  they  are  still  negative  for 
this  value  of  77.  Comparison  of  Figures  10  and  12  shows  that  mode  II  stress  intensity 
factors  decrease  as  77  decreases.  It  is  also  seen  that,  the  effect  of  the  nonhomogeneity 
parameter  on  the  stress  intensity  factors  is  quite  significant.  Again,  the  results  for 
7 d  =  0.0001  are  in  exact  agreement  with  the  results  obtained  from  the  homogeneous 
formulation.  In  Figures  13-14  and  15-16  modes  I  and  II  stress  intensity  factors  are  given 
for  77  =  0.4  and  77  =  0.8,  respectively.  It  can  be  observed  that  gradually  kj  becomes 
positive  and  kjj  becomes  negative  as  the  the  tangential  force  increases.  The  contact  stress 
distribution  is  shown  in  Figures  17-20  for  a  =  0.  In  this  case  the  stress  singularities  a  and 
j3  at  the  end  points  a  =  0  and  6  are  given  by  (56)  and  (174d),  respectively.  Figure  17 
shows  that  for  77  =  0  there  is  no  singularity  at  the  trailing  end  a  =  0.  In  Figures  18  and 
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19,  contact  stress  distributions  for  77  =  0.2  and  77  =  0.4  are  given.  For  these  relatively 
small  values  of  77  the  stress  singularities  at  b  are  greater  than  those  at  a  —  0  (i.e., 
—  >  —  a),  hence  the  skewed  distribution  in  Figures  18  and  19.  On  the  other  hand  for 

relatively  large  values  of  77,  |a|  >  \(5\  the  trend  is  reversed  and  there  is  a  greater  stress 
concentration  near  the  end  a  =  0  (see  Figure  20). 

Figures  21-28  show  the  modes  I  and  II  stress  intensity  factors  as  functions  of  a/ d  for 
a  constant  relative  contact  area  length  (b  —  a)/d  —  0.1.  These  figures  also  show  that,  the 
limiting  cases  of  7 d  —  0.0001,  are  in  very  good  agreement  with  the  results  obtained  from 
the  solution  of  the  homogeneous  half  plane  problem.  As  seen  in  Figure  21,  for 
7 d  =  0.0001  mode  I  stress  intensity  factors  are  negative  for  all  values  of  aj d,  which 
would  lead  to  crack  closure.  It  can  also  be  seen  that  mode  I  stress  intensity  factors  in  a 
graded  medium  are  larger  than  those  for  the  homogenous  medium  and  for  some  values  of 
7 d  and  a/d,  mode  I  stress  intensity  factors  become  positive.  Figure  22  shows  that  mode  II 
stress  intensity  factors  are  positive  for  all  values  of  a/d  in  a  homogeneous  medium,  and 
they  decrease  gradually  as  the  nonhomogeneity  parameter  7 d  increases.  Figures  23  and  24 
show  the  results  for  77  =  0.2.  As  the  coefficient  of  friction,  hence  the  tangential  force 
increases  mode  I  stress  intensity  factors  increase  and  mode  II  stress  intensity  factors 
decrease.  The  results  for  77  =  0.4  and  77  =  0.8  are  shown  in  Figures  25-26  and  27-28 
respectively.  Figure  27  shows  that  for  77  =  0.8  mode  I  stress  intensity  factors  are  positive 
for  all  values  of  7 d. 

Contact  stress  distribution  for  (b  -  a)/d  =  0.1  and  a/d  -  0.4  are  given  in  Figures 
29-32.  Figure  29  shows  the  results  for  77  =  0.  Although  there  is  no  tangential  force  and 
singularities  are  equal  at  both  ends  of  the  contact  area,  the  stress  distribution  is  not 
exactly  symmetric  due  to  the  effect  of  the  surface  crack  in  the  graded  medium.  It  can  be 
seen  in  Figures  30-32  that,  as  the  coefficient  of  friction  increases,  singularity  at  the 
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leading  end  i.e.,  -  (3  decreases  and  there  is  a  higher  stress  intensification  at  the  trailing 
end. 

Another  set  of  results  for  the  stress  intensity  factors  are  given  in  Figures  33-40  for  a 
relatively  larger  stamp  size  (b  -a)/d  =  1.0.  The  trends  are  similar  as  in  Figures  21-28. 
The  contact  stress  distributions  for  ( b  —  a)/d  =  1.0  are  shown  in  Figures  41-44  for 
various  values  of  the  friction  coefficient  rj. 

Some  conclusions 

1.  Analytically  the  contact  problem  for  a  graded  half-plane  with  exponentially  decaying 
stiffness  is  not  a  well-posed  problem. 

2.  The  trailing  end  of  the  sliding  rigid  stamp  with  friction  is  a  likely  location  of  surface 
crack  initiation  due  to  greater  stress  concentration. 

3.  In  the  medium  containing  a  surface  crack  and  loaded  by  a  sliding  rigid  stamp,  the 
mixed  mode  stress  state  at  the  crack  tip  is  such  that  the  cracks  tend  to  be  periodic  and 
curve  backward. 

4.  In  the  coupled  crack/contact  problems  for  a  graded  medium  stress  singularities  a,  (3 
and  to  are  independent  of  the  material  nonhomogeneity  constants  7  and  /io  and  depend  on 
the  friction  coefficient  rj  and  the  surface  value  of  the  Poisson's  ratio  (through  the  elastic 
constant  k)  only. 
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Figure  3:  The  distribution  of  the  contact  stress  on  the  surface  of  the  graded  medium 
loaded  by  a  flat  stamp  as  shown  in  Figure  2,  d  =  0,  rj  —  0,  k  =  2. 
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Figure  4:  The  distribution  of  the  in-plane  stress  on  the  surface  of  the  graded  medium 
loaded  by  a  flat  stamp  as  shown  in  Figure  2,  d  =  0,  r)  —  0,  k  =  2. 
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Figure  7:  The  distribution  of  the  contact  stress  on  the  surface  of  the  graded  medium 
loaded  by  a  flat  stamp  as  shown  in  Figure  2  ,  d  =  0,  rj  =  0.8,  k  =  2. 
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Figure  8:  The  distribution  of  the  in-plane  stress  on  the  surface  of  the  graded  medium 
loaded  by  a  flat  stamp  as  shown  in  Figure  2  ,  d  =  0,  rj  =  0.8,  k  =  2. 
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Figure  9:  Mode  I  stress  intensity  factors  for  an  edge  crack  in  a  graded  half  plane  indented 
by  a  flat  stamp  as  shown  in  Figure  2,  a  =  0,  rj  =  0,  v  =  0.25. 
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Figure  10:  Mode  II  stress  intensity  factors  for  an  edge  crack  in  a  graded  half  plane 
indented  by  a  flat  stamp  as  shown  in  Figure  2,  a  =  0,  77  =  0,  v  =  0.25. 


Figure  11:  Mode  I  stress  intensity  factors  for  an  edge  crack  in  a  graded  half  plane 
indented  by  a  flat  stamp  as  shown  in  Figure  2,  a  =  0,  rj  =  0.2,  u  =  0.25. 
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Figure  12:  Mode  II  stress  intensity  factors  for  an  edge  crack  in  a  graded  half  plane 
indented  by  a  flat  stamp  as  shown  in  Figure  2,  a  =  0,  77  =  0.2,  u  =  0.25. 
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Figure  13:  Mode  I  stress  intensity  factors  for  an  edge  crack  in  a  graded  half  plane 
indented  by  a  flat  stamp  as  shown  in  Figure  2,  o  =  0,  r]  =  0.4,  v  =  0.25. 


Figure  14:  Mode  13  stress  intensity  factors  for  an  edge  crack  in  a  graded  half  plane 
indented  by  a  flat  stamp  as  shown  in  Figure  2,  a  =  0,  rj  =  0.4,  u  =  0.25. 


Figure  15:  Mode  I  stress  intensity  factors  for  an  edge  crack  in  a  graded  half  plane 
indented  by  a  flat  stamp  as  shown  in  Figure  2,  a  =  0,  rj  =  0.8,  v  =  0.25. 


Figure  16:  Mode  II  stress  intensity  factors  for  an  edge  crack  in  a  graded  half  plane 
indented  by  a  flat  stamp  as  shown  in  Figure  2,  a  =  0,  rj  =  0.8,  u  —  0.25. 
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Figure  17:  Contact  stress  distribution  for  a  graded  half  plane  with  an  edge  crack  and 
indented  by  a  flat  stamp  as  shown  in  Figure  2,  a  =  0,  r]  =  0,  u  =  0.25,  b/d  —  0.5. 


-1.0 
0.0  i — ' 


^(o,  y) 

p/b  -i-0 


yd  =  0.0001 
yd  =  0.5 
yd=\.0 
yd  =  2.0 


Figure  18:  Contact  stress  distribution  for  a  graded  half  plane  with  an  edge  crack  and 
indented  by  a  flat  stamp  as  shown  in  Figure  2,  a  =  0,  rj  =  0.2,  v  =  0.25,  b/d  =  0.5. 
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Figure  19:  Contact  stress  distribution  for  a  graded  half  plane  with  an  edge  crack  and 
indented  by  a  flat  stamp  as  shown  in  Figure  2,  a  =  0,  r?  =  0.4,  v  =  0.25,  b/d  =  0.4. 
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Figure  20:  Contact  stress  distribution  for  a  graded  half  plane  with  an  edge  crack  and 
indented  by  a  flat  stamp  as  shown  in  Figure  2,  a  =  0,  rj  =  0.8,  v  =  0.25,  b/d  —  0.4. 
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Figure  21:  Mode  I  stress  intensity  factors  for  an  edge  crack  in  a  graded  half  plane 
indented  by  a  flat  stamp  as  shown  in  Figure  2,  (b  —  a)/d  =  0.1,  rj  —  0,  v  =  0.25. 


Figure  22:  Mode  II  stress  intensity  factors  for  an  edge  crack  in  a  graded  half  plane 
indented  by  a  flat  stamp  as  shown  in  Figure  2,  (b  —  a)/d  =  0.1,  rj  =  0,  v  —  0.25. 


Figure  23:  Mode  I  stress  intensity  factors  for  an  edge  crack  in  a  graded  half  plane 
indented  by  a  flat  stamp  as  shown  in  Figure  2,  (b  —  a)/d  =  0.1,  r)  =  0.2,  v  =  0.25. 


Figure  24:  Mode  H  stress  intensity  factors  for  an  edge  crack  in  a  graded  half  plane 
indented  by  a  flat  stamp  as  shown  in  Figure  2,  (b  —  a)/d  =  0.1,  17  =  0.2,  u  =  0.25. 
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Figure  25:  Mode  I  stress  intensity  factors  for  an  edge  crack  in  a  graded  half  plane 
indented  by  a  flat  stamp  as  shown  in  Figure  2,  ( b  —  a)/ d  =  0.1,  rj  =  0.4,  v  =  0.25. 


Figure  26:  Mode  II  stress  intensity  factors  for  an  edge  crack  in  a  graded  half  plane 
indented  by  a  flat  stamp  as  shown  in  Figure  2,  (6  —  a)/d  =  0.1,  rj  =  0.4,  v  =  0.25. 
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Figure  27:  Mode  I  stress  intensity  factors  for  an  edge  crack  in  a  graded  half  plane 
indented  by  a  flat  stamp  as  shown  in  Figure  2,  (b  -  a)/d  =  0.1,  rj  =  0.8,  u  =  0.25. 
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Figure  28:  Mode  II  stress  intensity  factors  for  an  edge  crack  in  a  graded  half  plane 
indented  by  a  flat  stamp  as  shown  in  Figure  2,  ( b  —  a)/d  =  0.1,  rj  =  0.8,  v  =  0.25. 
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Figure  29:  Contact  stress  distribution  for  a  graded  half  plane  with  an  edge  crack  and 
indented  by  a  flat  stamp  as  shown  in  Figure  2,  (b  —  a)/d  =  0.1,  r)  =  0,  u  =  0.25, 
a/d  =  0.4. 
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Figure  30:  Contact  stress  distribution  for  a  graded  half  plane  with  an  edge  crack  and 
indented  by  a  flat  stamp  as  shown  in  Figure  2,  ( b  —  a)/d  =  0.1,  rj  =  0.2,  v  =  0.25, 
a/d  =  0.4. 
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Figure  31:  Contact  stress  distribution  for  a  graded  half  plane  with  an  edge  crack  and 
indented  by  a  flat  stamp  as  shown  in  Figure  2,  (b  —  a)/d  =  0.1,  rj  =  0.4,  v  =  0.25, 
o/d  =  0.4. 
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Figure  32:  Contact  stress  distribution  for  a  graded  half  plane  with  an  edge  crack  and 
indented  by  a  flat  stamp  as  shown  in  Figure  2,  (b  —  a)/d  =  0.1,  rj  =  0.8,  v  =  0.25, 
a/d  =  0.4. 
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Figure  33:  Mode  I  stress  intensity  factors  for  an  edge  crack  in  a  graded  half  plane 
indented  by  a  flat  stamp  as  shown  in  Figure  2,  (6  —  a)/d  =  1.0,  r]  —  0,  v  —  0.25. 
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Figure  34:  Mode  II  stress  intensity  factors  for  an  edge  crack  in  a  graded  half  plane 
indented  by  a  flat  stamp  as  shown  in  Figure  2,  (b  —  a)/d  =  1.0,  rj  =  0,  u  =  0.25. 


Figure  36:  Mode  II  stress  intensity  factors  for  an  edge  crack  in  a  graded  half  plane 
indented  by  a  flat  stamp  as  shown  in  Figure  2,  (6  -  a)/d  =  1.0,  77  =  0.2,  u  =  0.25. 
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Figure  37:  Mode  I  stress  intensity  factors  for  an  edge  crack  in  a  graded  half  plane 
indented  by  a  flat  stamp  as  shown  in  Figure  2,  ( b  -  a)/d  —  1.0,  77  =  0.4,  v  =  0.25. 


Figure  38:  Mode  II  stress  intensity  factors  for  an  edge  crack  in  a  graded  half  plane 
indented  by  a  flat  stamp  as  shown  in  Figure  2,  ( b  -  a)/ d  =  1.0,  7?  =  0.4,  v  =  0.25. 
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Figure  39:  Mode  I  stress  intensity  factors  for  an  edge  crack  in  a  graded  half  plane 
indented  by  a  flat  stamp  as  shown  in  Figure  2,  ( b  —  a)/ d  =  1.0,  rj  =  0.8,  v  =  0.25. 
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Figure  40:  Mode  II  stress  intensity  factors  for  an  edge  crack  in  a  graded  half  plane 
indented  by  a  flat  stamp  as  shown  in  Figure  2,  (6  -  a)/ d  =  1.0,  77  =  0.8,  u  =  0.25. 
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Figure  41:  Contact  stress  distribution  for  a  graded  half  plane  with  an  edge  crack  and 
indented  by  a  flat  stamp  as  shown  in  Figure  2,  ( b  -  a)/d  =  1.0,  rj  =  0,  u  =  0.25, 
o/d  =  0.4. 
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Figure  42:  Contact  stress  distribution  for  a  graded  half  plane  with  an  edge  crack  and 
indented  by  a  flat  stamp  as  shown  in  Figure  2,  (b  -  a)/d  =  1.0,  r]  =  0.2,  v  =  0.25, 
a/d  =  0.4. 


(2  y-(b  +  a))/(b-a) 


Figure  43:  Contact  stress  distribution  for  a  graded  half  plane  with  an  edge  crack  and 
indented  by  a  flat  stamp  as  shown  in  Figure  2,  ( b  —  a)/d  =  1.0,  77  =  0.4,  v  =  0.25, 
a/d  =  0.4. 
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Figure  44:  Contact  stress  distribution  for  a  graded  half  plane  with  an  edge  crack  and 
indented  by  a  flat  stamp  as  shown  in  Figure  2,  (6  —  a)/d  =  1.0,  rj  =  0.8,  v  =  0.25, 
a/d  =  0.4. 
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APPENDIX  A 


Some  Useful  Integrals 

Here,  we  give  the  expressions  that  are  used  to  evaluate  the  integrals  involving  the 
asymptotic  expansions  of  the  integrands  of  the  kernels.  There  are  three  types  of  integrals 
involving  the  asymptotic  expansions.  The  expressions  for  each  type  are  given  below. 

Integrals  of  type  1 


In  this  case,  the  integrals  that  we  want  to  evaluate  are  in  the  following  form, 


f 00  l 

Cn  =  J  —cos  (pu)dp, 
f°°  l 

Sn  =  /  —sin  (pu)dp, 
Ja  P 


n  =  1,2,3.. 

.,N, 

(Fla) 

n  =  1,2, 3.. 

.,N. 

(Fib) 

For  n  =  1,  following  results  are  obtained  using  MAPLE, 


Cl  =  -  Ci(AM), 

Si  =sign(u)(|-Si(A|u|)), 


(F2a) 

(F2b) 


where  Ci()  and  Si()  are  cosine  and  sine  integrals,  respectively,  and  they  are  defined  by 


cos(a)  —  1 


.  /  cos(aJ 

Ci(x)  —  70  +  ln(a:)  +  /  - 

Jo  a 


da. 


Si(x)  =  r^a) 
Jo  O' 


da. 


(F3a) 

(F3b) 


and  the  Euler  constant  is  70  =  0.5772156649.  For  n  >  1,  integrating  (Fla)  and  (Fib)  by 
parts  the  following  general  recursive  relations  are  obtained: 
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C,  =  -  ^^  +  r^-,,  »>1 

1  —  n  An  1  i  —  n 


1  —  n 


Sn  — 


1  sin(uA)  u 
1  —  n  An_1  1  —  n 


n  >  1. 


Following  result  is  also  used  in  the  integration  of  asymptotic  expressions, 

f°°  sin(pw)  7 r  . 

/  dp  =  -sign(u). 

Jo  P  2 


Integrals  of  type  2 


In  this  case,  we  consider  the  following  integral: 


(F4a) 


(F4b) 


/•oo  1 

fln  =  /  —  exp(pu)dp, 

Ja  p 


=  1,2,3...,  JV,  [u  <  0 


For  n  =  1,  following  result  is  obtained  using  MAPLE, 


J?!  =  Ei(  -  Au),  Ei(2:)  =  -  Ei(  -  z), 


(F7a,b) 


where  Ei()  is  the  exponential  integral  function.  For  n  >  1,  following  expression  is  used 
which  is  given  by  Gradshteyn  and  Ryzhik  [10]: 

f°°  exp(  -  px)dx  .  n+1  p^Ei(  — pA)  exp(  -  pA)  (  -  1  )fepfcAfc 


=  (-l)r 


—  1)...  (n  —  k)’ 


[P>  0 


(F8)  reduces  to  (F7)  for  n  =  1. 


Integrals  of  type  3 


Type  3  integrals  are  in  the  following  form: 
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f°°  1 

Cn  =  Ja  —cos(pv)exp(pu)dp, 


Sn  = 


-sm(pv)exp(pu)dp, 


=  1,2, 3...  ,iV, 

[n<0]. 

(F9a) 

=  1, 2,3...,  JV, 

[«<0]. 

(F9b) 

(F9)  are  evaluated  using  the  following  expressions  which  are  given  by  Gradshteyn  and 
Ryzhik  [10]: 


*exp(  —  /3x)cos(6x)dx  —  -(/?  4-  iS)  *T (p,((3  +  i6)A) 

£ 

+  \  (0  -  isrrifi,  (0  -  iS)A),  [8(0)  >  |3 (6)  |] : 


(F10a) 


1exp(  —  /3x)sin(Sx)dx  =  -(/?  +  i6)  ^(p,  (/5  +  i6)A) 

A 

-  i(0-«)-'TO*,(0-tf)A),  [8(0)  >  |9(f)|], 


(FlOb) 


where  T(, )  is  the  incomplete  Gamma  function.  Following  result  is  also  used  in  the 
integration  of  asymptotic  expressions, 


f°°  i 

/  — sin(pu)exp(pu)dp  =  —  arctan 

Jo  P 


u  <  0. 


(FI  1) 
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APPENDIX  B 


Some  leading  terms  in  asymptotic  expansions 

Here,  we  give  the  leading  terms  of  asymptotic  expansions  K$f°(u,y),  K^(u,y)  and 
y)i  y )  which  are  given  by  equations  (129)  and  (147),  respectively. 

cos(7y/2)(exp(  -  %)  -  l)exp(«$iy/2) 


g*n=-2 
9n  =  2 


*ii=  -2 


1* 

‘21 


(1  +  fv)<5l 

sin(7i//2)(exp(  -  fryy)  -  l)exp(%/2) 

(1  +  k)6\ 

sin(7y/2)(exp(  -  fry)  -  l)exp(fty/2) 
(1  +  k)6i 

cos(7y/2)(exp(  -  6yy)  -  l)exp(^y/2) 

(1  +  k)6  i 


-2 


(Gla) 

(Gib) 

(Glc) 

(Gld) 
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APPENDIX  C 


Kernels  of  the  integral  equations 


In  this  Appendix,  we  give  the  transformed  form  of  the  kernels  that  are  used  in  equations 
(182a,c)  and  the  terms  that  are  used  in  equations  (185). 


Hn{si,r)  =  ^{HUs{si,r)  +  Hnf(si,r)),  (Hla) 

Hn(si,r)  =  +  (Hlb) 

H22(s2,r)  =  ^(H22s{s2,r)  +  H22f(s2,r)),  (Hlc) 

H23(si,r)  =  — — — (H23s(s2,r)  +  Hi3f(s2,r)),  (Hid) 

Hzi(s3,r)  =  ^{H3u{s3,r)  +  Hzif(s3,r)),  (Hie) 

H32(s3,r )  =  ^(H32s{s3,r)  +  Hz2f(s3,r)),  (Hlf) 

Hzz(sz,r)  =  b-^Hzzf(sz,r),  (Hlg) 

where, 


Hijs^Si)  r)  —  /ljj5(x,t), 


(H2a) 


Hijf(shr)  —  hijf(x,t), 


(  d  d 

2Si+r 


x  —  \ 


b  —  CL  b  +  CL 

—  5i+— , 


i  =  1,2 


i  =  3 


(H2b) 


(H3a) 
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.7  =  1,2 


r  d  d 

2r  +  2’ 


t  =  < 


b  —  a  b  +  a  . 

.r+__)  j  =  3 


2  2 

The  terms  used  in  equations  (185)  are  given  below: 

2ai-1/2r(  -  l/2)T(n  +  +  1) 

= - irVin  +  1/2 +  ai) - X 

x  F(n+1,  -  n  +  1/2  -  on; 3/2;  (1  -  s^/2) 

+  j\  1  -  r)-1/2(l  +  r)ipW/3«>(r)ffn(»,,r)* 
m13„(Sl)  =  /  (1  -  r)'(l  +  r)“!pW“!>(r)if13(Sl,r)dr, 


(H3b) 


(H4a) 

(H4b) 


^22n(S2)  = 


2ai~1/2r(  -  l/2)r(n  +  Qi  +  1) 

7iT(n  +  1/2  +  cci) 


x  f(tz+1,  -  ra  +  1/2  -  <*i;  3/2;  (1  -  s2)/2)  + 


+/ 

(1  -  r)-1/2(l  +  r)“*P<-122'“'>(r)P11(s2,  r)dr, 

-1 

(H4c) 

^23n(S2)  = 

yV  -  r/(l  +  r)“2Pf  ’ Q2)(r)^13(s2,  r)dr, 

(H4d) 

m3ln(S3)  = 

^(1  -  r)-1/2(  1  +  rpPt1/2'ai)(r)H31(S!i, r)dr, 

(H4e) 

^32^(53)  = 

J  (1  -  r)_1/2(l  +  r)“1P7^_1/2’Q:i)(r)iT3i(s3,r)dr, 

(H4f) 
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